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ABSTRACT 

Loran-C  time  difference  (TD)  readings  were  taken 
consecutively  at  several  survey  stations  near  Monterey,  CA, 
over  a  period  of  three  days.  One  station,  Range-7,  was 
designated  the  "known"  point  of  a  differential  Loran-C 
system.  Readings  from  the  known  point  were  used  to  correct 
readings  from  a  second  survey  station,  Luces  Point.  A  method 
of  improving  the  precision  of  Loran-C  TD  readings  based  on 
the  redundancy  and  relative  accuracy  of  three  LOP ' s  was 
developed  and  applied  to  the  data.  Since  only  one  receiver 
was  available,  a  linear  regression  of  TD  vs  time  was 
calculated  and  used  in  the  differential  correction. 

Based  on  496  sets  of  data  taken  at  5-second  intervals  at 
Range-7  in  two  groups,  before  and  after  about  250  readings  at 
Luces  Point,  the  absolute  accuracy  of  Luces  Point  data  was 
improved  from  about  385  m  to  about  48  m  compared  to  the  known 
position  of  the  point.  Precision  was  improved  from  about 
14.9  m  to  about  12.6  m  circular  error  (CEP)  using  the  three- 
station  correction.  Further  improvement  would  probably  have 
resulted  if  two  receivers  were  available  for  real  time 
corrections . 
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I.      INIRXUCTCCN 

A.     Purpose 

The  purpose  of  this  study  was  to  investigate  the  potential 
accuracy  of  the  Loran-C  navigation  system  using  a  differential 
correction  technique  and  a  correction  for  the  use  of  three  secondary 
stations  rather  than  the  usual  two. 

Sources  of  error  were  first  examined  and  their  magnitudes 
estimated.  Methods  of  expressing  the  accuracy  of  a  determination  of 
geographic  position  were  also  considered. 

It  was  determined  that  position  accuracy  is  composed  of  two 
factors:  the  absolute  accuracy  of  the  position  (which  can  be 
considered  the  bias,  or  systematic  error,  in  the  system)  and  the 
precision  of  the  individual  readings. 

Over  longer  periods  of  data  collection  it  was  considered  possible 
to  subtract  out  the  system  bias  by  means  of  a  quasi-differential 
correction  based  on  intermittent  time  difference  readings  at  a  fixed 
point,  designated  the  "known"  point  of  the  differential  system,  and 
to  improve  the  precision  of  TD  readings  by  taking  advantage  of  the 
redundancy  of  the  reception  of  three  master-secondary  pairs. 

For  this  purpose,  a  theory  of  corrections  was  developed  and  an 
experiment  was  designed  to  determine  the  effectiveness  of  the 
corrections.  The  experimental  data  was  processed  by  computer,  and  the 
relative  accuracies  were  calculated.  Computer  programs  for  processing 
and  displaying  data  were  written  for  this  purpose. 
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B.  Method 

To  collect  data  for  processing,  several  known  survey  points  with 
published  positions  were  chosen  as  locations  for  collecting  Loran-C 
data.  Using  known  points  provided  a  means  of  comparing  results.  One 
point  in  a  relatively  central  location  was  chosen  as  the  "known"  point, 
where  a  larger  amount  of  data  was  collected.  Data  from  the  other  points 
was  corrected  based  on  the  data  from  the  "known"  point. 

During  periods  of  about  eight  hours,  the  receiver  was  moved  from 
one  survey  point  to  another.  The  antenna  was  set  up  over  the  points 
and  time  difference  data  was  collected  for  periods  of  about  20  to  30 
minutes  at  each  point.  The  time  required  to  move  the  receiver  from  one 
point  to  another  varied  from  about  30  to  60  minutes. 

Data  was  collected  by  a  Racal  Megapulse  monitor  receiver  with 
internal  clock  accuracy  of  about  0.02  ^isec.  A  Silent  700  data  terminal 
was  used  to  collect  the  data  on  cassette  tapes.  Data  was  recorded  at 
intervals  of  5  seconds.  Later,  the  data  was  transferred  to  5V  disks 
for  use  on  a  computer. 

Since  there  was  no  way  to  collect  time  difference  data  at  two 
points  simultaneously,  a  method  of  estimating  the  time  differences  at 
the  known  point  was  developed.  This  involved  using  the  available  data 
to  establish  linear  regression  lines  of  time  difference  versus  time  at 
the  known  point. 

Using  the  known  time  of  data  collection  at  the  other  points,  an 
estimate  of  the  time  differences  at  the  known  point  could  be  used  to 
calculate  a  correction  for  data  at  the  other  points  at  the  times  of 
data  collection. 
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A  number  of  computer  programs  were  written  to  calculate  the 
required  information  for  differential  corrections.  A  Correction  For 
Three  Secondary  Stations  was  developed  mathematically  and  used  to 
process  the  data  in  an  effort  to  reduce  the  random  error  of  the  time 
difference  readings. 

Based  on  reduction  of  standard  deviations  of  the  data  and 
reduced  circular  error,  the  three  station  correction  was 
effective  to  some  degree,  although  absolute  error  turned 
out  to  be  the  limiting  factor  in  the  experiment.  Reduction 
of  absolute  error  by  means  of  the  differential  correction 
was    significant. 


II.   THE  LCRAN-C  NAVIGATION  SYSTEM 

A.  Basic  Description 

Loran-C  is  an  electronic  positioning  system  based  on  the  reception 
of  pulsed  LF  (90-110  kHz)  signals  from  a  master  transmitting  station 
and  one  or  more  secondary  stations  which  comprise  a  Loran-C  "chain".  A 
brief  description  of  the  system  will  be  provided  here  as  a  basis  for 
later  discussions  of  the  accuracy  considerations  and  differential 
methods  used  in  the  experiment. 

In  the  hyperbolic  mode,  geographic  position  (GP)  is  determined  from 
time  difference  (TD)  readings  provided  by  a  Loran-C  receiver.  The 
difference  in  reception  time  of  pulses  received  from  the  master  station 
and  pulses  received  from  each  secondary  station  give  the  receiver's 
position  in  terms  of  a  hyperbolic  lattice  system.  Geographic  position 
can  be  calculated  from  the  TD  readings  using  Loran-C  lattice  tables 
published  by  the  Defense  Mapping  Agency  or  computed  directly  using 
various  computer  techniques.  Nautical  charts  overprinted  with  Loran-C 
TD  lattices  are  also  available  for  plotting  geographic  position 
directly  without  converting  first  to  latitude  and  longitude. 

The  hyperbolic  lattice  system  is  a  series  of  hyperbolas.  A  pulse 
(actually  a  group  of  pulses)  is  first  transmitted  by  the  master 
station.  The  wave  travels  to  each  secondary  station  of  the  Loran-C 
chain,  and  each  secondary  transmits  its  own  pulse  after  a  time  delay 
(called  the  coding  delay)  which  is  unique  to  each  secondary. 


The  master  pulse  is  received  first  by  the  receiver.  A  short  time 
later  the  pulses  from  the  secondaries  are  received.  The  tijme  dif- 
ference between  reception  of  the  master  and  each  secondary  is  measured 
in  the  receiver.  Each  TD  identifies  one  line  of  position  (LDP)  in  the 
set  of  IDP's  for  a  master-secondary  pair.  Figure  II-l  shows  a  typical 
Icran-C  master-secondary  pair  with  a  set  of  hyperbolic  LDP's.  Each  LDP 
in  this  set  represents  a  line  of  equal  time  difference  between 
reception  of  the  master  pulse  and  reception  of  the  pulse  from  the 
secondary.  Figure  II-2  shows  an  example  of  a  hyperbolic  lattice  system 
formed  by  two  master-secondary  pairs. 

As  with  latitude  and  longitude,  two  LDP's  are  required  to  define  a 
geographic  position.  Unlike  latitude  and  longitude,  two  TD  readings 
from  different  secondaries  do  not  uniquely  define  a  GP,  since  the  TD's 
may  intersect  at  two  points  within  the  coverage  area.  This  is  normally 
not  a  problem,  since  even  an  approximate  knowledge  of  the  geographic 
position  of  the  receiver  will  allow  determination  of  the  correct 
intersection  point. 

When  more  than  two  LDP's  are  available  to  define  a  GP  (when  signals 
from  more  than  two  secondary  stations  can  be  received) ,  the 
intersection  of  the  three  LDP's  define  three  geographic  points  (see 
Figure  LV-2)  .  The  accuracy  of  the  two  best  pairs  of  IDP's  is 
sufficient  for  most  navigation  purposes,  but  it  is  possible  to  combine 
three  or  more  LDP's  to  estimate  the  position  of  one  point.  A  method 
using  corrections  from  three  LDP's  is  used  in  the  experiment  and  will 
be  discussed  in  Section  IV. B. 


FIGURE  II-l. 


Typical  set  of  hyperbolic  LDP's  associated  with  a  master- 
secondary  pair. 


FIGURE  II-2, 


Typical  hyperbolic  lattice  formed  by  sets  of  LDP's 
associated  with  two  master-secondary  pairs  with  a 
common  master  station. 


B.  Time  Dif ference  Equations 

The  difference  between  arrival  time  of  the  signal  from  the  master 
and  arrival  time  of  the  signal  from  a  secondary  can  be  described  by 
the  following  equations  (Ref .  1) : 

The  basic  time  difference  equation  for  Loran-C  is: 

TD  =  tg  -  tn  [II-l] 

where: 

TD  =  time  difference  at  a  given  point  P. 
tg  =  time  of  reception  of  secondary  pulse, 
tm  =  time  of  reception  of  master  pulse. 

A  more  useful  form  of  this  equation  takes  into  account  the 
components  which  make  up  the  travel  time  of  the  signals  over  each 
section  of  the  paths  from  the  stations  to  point  P: 

TD  =   (Ts  +  DTS  +  ED)  -  (Tm  +  DTm)  +  ASF         [II-2] 

where: 

TD  =  Time  Difference  at  a  given  point  P. 

Ts  =  Distance  in  time  units  between  the  secondary  station 
and  point  P. 

Tm  =  Distance  in  time  units  between  the  master  station 
and  point  P. 

DTS  =  Secondary  Phase  Correction  for  all  seawater  path 
between  secondary  station  and  point  P. 

DTm  =  Secondary  Phase  Correction  for  all  seawater  path 
between  master  station  and  point  P. 

ED  =  Emission  Delay  for  the  secondary  station. 

ASF  =  Additional  Secondary  Phase  Correction  for  point  P. 


Equation  II-2  contains  two  corrections  which  are  applied  to  Loran-C 
TD  readings.  The  velocity  of  electromagnetic  waves  in  air  under 
standard  conditions  is  changed  when  the  wave  propagates  over  a  surface. 
The  Secondary  Phase  Correction  Factor  (DTS  and  DIm  in  equation  II-2)  is 
a  time  correction  for  the  change  in  velocity  caused  by  propagation  of 
electromagnetic  waves  over  seawater.  This  change  in  travel  time  can  be 
calculated  by  the  following  empirically  derived  equations  (Ref .  2) : 

For  distances  greater  than  100  statute  miles  (T  >  537  /xsec) : 

DT  =  Aq/T  +  A1  +   A2*T  [H-3] 

where: 

DT  =  Secondary  Phase  Correction  for  all  seawater  path. 

Aq  =  129.04398  /xsec2 

A!  =  -0.40758  Msec 

A2     =  0.00064576438 

T       =  Travel  time   (Ts  or  Tm  in  equation  II-l)    in  /xsec. 

For  distances  less  than  100  statute  miles   (T  <  537  /xsec) : 

DT     =     Bq/T  +  Bx  +  B2*T  [H-4] 

where: 

Bo     =     2.7412979  usee2 

B;l     =     -0.011402  /Lisec 

B2  =  0.00032774624 

T   =  Travel  time  in  usee. 


In  equation  II-2  the  Additional  Secondary  Phase  Factor  (ASF) 
corrects  for  propagation  over  land  (actually  any  other  non-seawater) 
surfaces,  which  have  a  wide  variety  of  conductivities  (see  Section 
III-C-2  or  Ref .  3  for  further  explanation  of  ASF  Corrections) .  The  ASF 
correction  is  in  addition  to  the  Secondary  Phase  Correction  for  all 
seawater  path. 

Secondary  Phase  Corrections  for  all  seawater  path  are  always  added 
in  before  calculation  of  the  geographic  position  from  TD's.  ASF 
corrections  are  much  smaller,  and  are  often  determined  from  tables 
which  require  a  dead  reckoning  position  within  10'  of  latitude  and 
longitude  of  the  final  position.  If  no  DR  position  is  available,  the 
position  calculated  from  TD's  corrected  only  for  all  seawater  path  is 
used  as  a  DR  position  to  determine  ASF  corrections.  When  used  in 
hyperbolic  mode,  the  ASF  correction  in  Equation  II-2  is  a  net 
correction  which  applies  to  the  particular  master-secondary  pair  (the 
difference  between  the  range-range  corrections  for  each  station) . 

Emission  Delay  (ED)  in  Equation  II-2  includes  a  coding  delay,  which 
is  the  time  between  reception  of  the  master  pulse  at  the  secondary 
station  and  transmission  of  the  secondary's  pulse,  plus  one  way  travel 
time  from  the  master  to  the  secondary,  including  the  Secondary  Phase 
Correction  for  all  seawater  path  for  the  master-secondary  baseline  path 
(Ref.  4) . 

The  result  of  application  of  these  corrections  to  the  basic  time 
difference  between  reception  of  master  and  secondary  signals  is  a 
reduction  of  the  time  difference  to  a  path  through  air  with  an  Index  of 
Refraction  of  1.000338  (Ref.4:p.4). 
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C.  Differential  Corrections 

Differential  techniques  can  be  used  to  improve  the  absolute 
accuracy  of  Loran-C  TD  readings.  Differential  Loran-C  uses  monitor 
receivers  located  at  fixed  points  in  or  near  the  area  where  improvement 
is  needed.  The  most  common  method  is  to  use  TD  readings  of  the 
monitor  receivers  to  determine  the  magnitude  of  the  variations  in  the 
TD's  from  the  long  term  average  of  the  TD's  at  the  fixed  points.  These 
variations  are  then  applied  to  receiver  readings  in  the  area  of 
interest  before  geographic  positions  are  calculated  from  them. 
Corrections  are  updated  periodically  as  explained  later  in  this 
section.  (Ref .l:p.434) 

The  differential  correction  (including  geographic  correction)  to 
Loran-C  readings  at  an  unknown  point  P2  with  a  monitor  receiver  at  a 
known  point  P^   can  be  described  by  the  following  equation: 

TD  at  P2  =  Received  TD  at  P2  +  Differential  Cor.  +  Geographic  Cor. 

or 
TD2c  =  TD2  +  (TDX  -  TD1C)  +  TDg  [II-5] 

where 

TE>2c  =  Corrected  TD  reading  at  Point  P2  (the  unknown  point) . 
TD1C  =  Long  term  average  TD  at  P^  (location  of  the  monitor) . 
TDj   =  TD  reading  of  receiver  at  the  known  point  (P-J  . 
TD2   =  TD  reading  of  the  receiver  at  the  unknown  point  (P2) . 
TDg   =  TD  variation  which  is  a  function  of  location. 
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For  practical  purposes  it  is  possible  to  use  any  correction  from 
real  time  to  a  correction  averaged  over  a  period  of  up  to  several  days. 
The  method  used  in  this  experiment  (described  in  Section  IV-C)  was  to 
determine  a  regression  line  of  monitor  receiver  TD  readings  at  the 
known  point.  A  computed  correction  was  applied  to  the  TD  reading  at 
the  unknown  point  based  on  the  time  of  the  TD  reading  at  the  unknown 
point  and  the  equation  of  the  regression  line.  This  method  was  used 
because  only  one  monitor  receiver  was  available  and  it  was  used  as  both 
a  monitor  at  the  known  point  and  at  the  unkncwn  points.  In  this 
experiment  a  correction  based  on  reception  of  three  secondary  stations 
(rather  than  two)  was  applied  to  the  TD  readings  from  each  secondary 
station  before  the  regression  line  parameters  were  calculated,  as 
explained  later  in  Section  IV-B. 

The  geographic  correction  -  that  part  of  the  difference  in  TD 
readings  at  points  P-^  and  P2  which  is  a  due  to  the  differences  in  ASF 
corrections  at  the  two  points  -  is  a  particular  problem.  For 
geographic  points  at  sea,  tables  of  ASF  corrections  within  10'  of 
latitude  and  longitude  are  available,  and  may  be  accurate  enough  to 
provide  corrections  within  the  necessary  accuracy.  Tables  of 
corrections  for  geographic  points  on  land  are  not  published,  but  were 
obtained  for  this  experiment  from  the  Defense  Mapping  Agency 
Hydrographic/Topographic  Center. 

Test  results  of  differential  Loran-C  on  the  East  Coast  U.S.  chain 
have  shown  accuracies  of  about  4.6  m  CEP  at  distances  up  to  140  km  from 
the  monitor  site  (Ref .l:p.434) . 
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IH.   GENERAL  ANALYSIS  OF  ERRORS 

To  quantify  the  amount  of  error  and  compare  the  accuracy  of 
different  types  of  corrections  made  to  Loran-C  geographic  position 
fixes,  it  is  necessary  to  provide  some  basic  definitions.  The 
following  sections  give  the  definitions  used  in  this  study  for  the 
types  of  accuracy  and  measures  of  accuracy.  One  and  two  dimensional 
sources  of  error  in  Loran-C  positioning  are  listed  and  approximate 
values  for  the  error  are  given  where  possible.  Error  reduction  based 
on  differential  correction  is  described. 

A.  Types  of  Accuracy 

There  are  several  ways  to  specify  the  accuracy  of  electronic 
navigation  systems.  The  following  definitions  of  three  types  of 
accuracy  have  been  found  useful  in  this  study  (from  Ref .  5) : 

Absolute  Accuracy;  For  purposes  of  this  study,  absolute  accuracy 
is  a  measure  of  the  ability  to  determine  geographic  position  (in  terms 
of  latitude  and  longitude  or  other  fixed  coordinate  system)  from  a 
given  set  of  Loran-C  coordinates  (Ref .5:p.l72) .  The  error  includes  all 
types  of  systematic  and  random  error  associated  with  the  system. 

Absolute  accuracy  in  this  experiment  refers  to  the  accuracy  of  the 
calculation  of  the  positions  of  the  ,lunknownl,  points  based  on  receiver 
TD  readings,  three  station  corrections,  and  differential  corrections, 
compared  to  the  known  published  positions  of  those  points. 
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Repeatable  Accuracy;  The  common  navigation  definition  of 
repeatable  accuracy  is  the  ability  to  return  to  the  same  point  using  a 
given  set  of  electronic  coordinates  (Ref  .5:p.l72)  .  Error  is  measured 
in  linear  units.  This  differs  from  absolute  accuracy  in  that  the 
position  itself  need  not  be  known  in  terms  of  some  other  coordinate 
system.  Repeatable  accuracy  is  normally  better  than  absolute  accuracy 
because  the  part  of  systematic  error  which  is  a  function  of  position 
has  been  removed  (the  position  remains  the  same) ,  leaving  only  random 
error  and  time  dependant  systematic  error.  Repeatable  accuracy  in  the 
context  of  this  experiment  is  a  measure  of  the  precision  of  measurement 
of  the  TD's,  which  is  most  affected  by  the  electronic  circuitry  of  the 
transmitters  and  receivers  involved  (over  short  periods  of  time  in 
which  terrain  and  atmospheric  factors  do  not  change  significantly) . 

Repeatable  accuracy  is  similar  to  the  accuracy  expected  from 
differential  corrections.  The  difference  between  the  two  is  that 
since  repeatable  accuracy  refers  to  the  ability  to  return  to  the  same 
point  using  the  same  set  of  received  TD  readings,  there  is  no 
difference  in  the  ASF  (terrain)  correction  but  there  may  be  some 
difference  due  to  variability  of  TD's  with  time. 

Differential  corrections  have  different  ASF  corrections 
at  different  points.  These  must  be  calculated  and  applied 
to  the  TD  readings.  At  a  single  point  the  ASF  corrections 
will  probably  not  change  over  the  time  period  of  interest, 
but  over  longer  periods  some  error  will  be  introduced.  The 
uncertainty  in  ASF  corrections  could  be  a  significant 
source  of  error. 
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Relational  Accuracy:  The  accuracy  with  which  it  is  possible  to 
determine  the  position  of  one  observer  in  relation  to  another  observer 
using  the  same  system  (Ref .5:p.l72) .  A  form  of  relational  accuracy 
enters  into  differential  corrections  when  two  different  receivers  are 
used.  Even  at  the  same  point  at  the  same  time  there  will  usually  be  a 
slight  difference  in  the  Loran-C  readings  of  different  receivers.  This 
is  caused  by  errors  associated  with  the  receiver  circuitry  and  random 
errors  over  the  slight  differences  in  the  time  at  which  the  readings 
are  made  by  each  receiver,  since  all  other  factors  are  the  same.  The 
same  type  of  error  exists  with  two  different  receivers  separated  by 
some  distance  (as  with  differential  corrections) ,  although  this  type  of 
error  is  usually  small  compared  to  other  sources  of  error. 

B.  Measures  of  Accuracy 

In  addition  to  the  types  of  accuracy  listed  above,  it  is  necessary 
to  specify  how  accuracy  is  to  be  measured.  This  can  be  viewed  as  an 
estimate  of  the  amount  of  error.  The  following  are  definitions  of  the 
basic  methods  used  here  of  specifying  the  measures  of  accuracy  of 
position.  Their  usefulness  and  limitations  in  this  application  are 
also  described. 

Root  Mean  Square  (RMS)  Error:  This  is  the  linear  (one  dimensional) 
error  in  the  location  of  an  LDP  (Ref  .5:p.l72) .  It  can  be  specified  in 
linear  units  or  time  units  (meters  and  microseconds  are  used  here) . 
The  normal  specification  of  RMS  error  is  the  number  of  meters  or 
microseconds  representing  some  whole  number  of  standard  deviations  from 
the  mean.   In  the  application  used  here,  the  variation  in  TD  readings 
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from  each  master-secondary  pair  can  be  expressed  as  a  standard 
deviation  in  microseconds,  which  corresponds  to  a  displacement  of  the 
LDP  by  a  corresponding  number  of  meters  measured  normal  to  that  line  of 
position. 

Circular  Error  Probable  (CEP) :  This  is  the  two  dimensional  error 
which  specifies  the  radius  of  a  circle  within  which  there  is  a  50% 
probability  of  being  located,  based  on  positions  calculated  from  the  TD 
readings  of  the  receiver.  The  estimated  position  of  the  point  (based 
on  latitude  and  longitude  or  other  coordinate  systems)  and  the  radius 
of  the  circle  are  used  to  specify  this  measure  of  accuracy.  Circles 
other  than  the  standard  50%  circle  can  be  used  (67%  and  90%  are 
common).  (Ref .5:p.l73)  The  main  advantage  of  using  CEP  (rather  than 
the  error  ellipse  described  below)  is  that  direct  comparison  between 
accuracies  of  different  systems,  techniques,  or  at  different  points  can 
be  made.  CEP  in  this  experiment  is  used  as  a  measure  of  the  precision 
of  the  TD's  indicated  by  the  receiver  at  a  fixed  point  and  the 
precision  of  corrected  TD's  based  on  these  uncorrected  TD's. 
Corrections  for  the  absolute  position  (the  bias  of  the  TD's  from  true 
readings)  is  of  secondary  importance  in  this  application. 

Circular  error  estimates  in  this  experiment  are  based 
on  the  radius  which  50%  of  the  data  are  greater  and  50%  are 
less,  as  calculated  by  a  sorting  program. 
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Error  Ellipses:  Error  ellipses  are  the  preferred  way  of 
describing  the  accuracy  of  electronic  navigation  systems  for  many 
purposes.  As  with  circular  error,  accuracy  is  specified  in  terms  of 
the  probability  of  an  estimate  of  geographic  position  falling  within  a 
geometric  figure  centered  on  the  actual  geographic  position  (an  ellipse 
rather  than  a  circle  is  used)  (Pef  .6) .  Unlike  the  circle  used  in  CEP, 
an  ellipse  provides  additional  information  about  the  directional 
characteristics  of  the  accuracy,  which  can  be  very  important  in  some 
applications.  The  specifications  for  elliptical  accuracy  are  usually 
the  lengths  of  the  major  and  minor  axes  and  the  orientation  of  the 
ellipse  (such  as  the  azimuth  of  the  direction  of  the  major  axis) .  A 
disadvantage  of  using  error  ellipses  in  this  experiment  is  the 
difficulty  of  directly  comparing  the  accuracy  of  different  ellipses 
with  different  orientations.  For  that  reason  most  of  the  accuracy 
comparisons  in  this  study  have  been  specified  in  terms  of  circular 
error.  Error  ellipses  have  been  used  where  appropriate. 

Error  ellipses  are  related  to  the  orientation  of  two  intersecting 
LDP's  by  the  following  equations  (Ref .6:p.72-74) : 


a^2  sin  20 


tan  29  =  [III-l] 

a!2 (cos  20)  +  a22 


°12  +  °22   +  ((al2  +  a22)2  "  4(sin20)a12a22)^ 
ax2  = [IH-2] 


2  sin20 
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°12   +  °22  ~    ((°12  +  °22)2   "  4(sin20)a12a22)"2 
av2  =  [III-3] 

2  sin20 


'y 


where: 

G      =      the    angle    of    the    major    axis    of    the    ellipse 
from   the    direction    of    the    LOP    with    smaller 
variance,    toward    the    LOP    with    greater    variance, 
in    the    direction    of    the    smaller    angle. 

o-^  =    variance  of  the  LOP  with  smaller  variance. 

o~22  =    variance  of  the  LOP  with  the  greater  variance. 

0      =      the    smaller    angle    between    LOP ' s . 

These  equations  were  used  to  calculate  the  dimensions  and 
orientation  of  the  error  ellipse  for  Range-7,  the  "known"  point  used 
for  differential  corrections    (see  Figure  V-3) . 

C.     Sources  of  Linear  Error  in  Loran-C  Time  Differences 

The  position  fixing  accuracy  of  the  LORAN-C  system  depends  on  a 
number     of     factors.  These     can     be     divided     into     two     categories 

(Ref .7:p.l8) :  temporal  errors  affecting  the  repeatability  of  the 
system  (errors  which  vary  with  time) ,  and  systematic  (time  invariant) 
errors.  Both  of  these  types  of  error  may  also  be  functions  of 
geographic  position. 

The  effect  of  these  errors  is  to  produce  TD  readings  at  the 
receiver  which  differ   from  those  which  would  be   expected   on  a  perfect 
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hyperbolic  lattice  on  the  surface  of  a  perfect  ellipsoid,  resulting  in 
error  in  the  calculation  of  the  geographic  position.  A  linear  error  of 
1.0  jLtsec  in  time  difference  is  approximately  equivalent  to  150  m  on  the 
baseline,  with  one  dimensional  error  at  other  points  being  up  to  about 
ten  times  greater,  depending  on  the  location  within  the  geometry  of  the 
system. 

The  main  sources  of  error  in  the  determination  of  TD's  used  to 
calculate  geographic  position  can  be  classified  as  follows  (Ref .8) : 

•  Transmitter  synchronization  timing  errors. 

•  Errors  associated  with  propagation  effects. 

•  Atmospheric  noise. 

•  Receiver  errors. 

1.  Transmitter  Synchronization  Timing  Errors 

Since  time  differences  between  reception  of  a  master  pulse  and 
secondary  pulses  are  used  to  calculate  geographic  position  in 
hyperbolic  Ioran-C,  it  is  essential  that  the  synchronization  of  timing 
between  transmissions  from  the  master  and  secondaries  be  highly 
precise . 

In  order  to  accomplish  this,  each  station  is  individually  timed 
using  cesium  clocks  which  drive  pulse  and  group  repetition  interval 
timing  circuits.  Independent  monitor  receivers  are  used  to  provide 
data  to  computers  which  adjust  the  relative  timing  between  master  and 
secondaries  of  the  chain  in  small  steps  (usually  10  to  20  ns) . 
Absolute  timing  of  the  master  stations  is  adjusted  to  operate  in 
coordination  with  Universal  Time,  Coordinated.  (Ref  .8 :p.  1129) 


18 


Estimates  of  timing  error  due  to  transmitter  synchronization 
are  on  the  order  of  0.05  /xsec  (Ref  .9:p.537) .  Actual  data  from  an 
operational  Loran-C  chain  showed  time  differences  controlled  within 
0.024  /isec  of  the  standard  value,  with  a  standard  deviation  (in 
absolute  time)  of  0.032  /usee  (Ref  .10:p.227) .  A  transmitter 
synchronization  error  of  0.01  usee  is  believed  to  be  attainable  if 
ground  monitor  data  is  used  to  remove  the  mean  synchronization  error 
(Ref.ll:p.401). 

2.  Errors  Associated  with  Propagation  Effects 

This  includes  most  types  of  error  external  to  the  electronics 
of  the  system,  and  is  the  main  source  of  error  in  the  determination  of 
geographic  position. 

Factors  which  influence  groundwave  propagation  include  the 
following  (Ref .12: pp. 39-53  and  Ref .13: pp. 173-187) : 


•  Electrical  Conductivity  of  the  surface  material 

•  Dielectric  Constant  of  the  surface 

•  Index  of  Refraction  of  the  air 

•  Lapse  Rate  of  the  Index  of  Refraction  with  altitude 


These  factors  are  complex  functions  of  several  variables,  the 
main  ones  being: 


•  type  and  composition  of  the  surface  over  which  the  wave  travels 

•  topography  of  the  surface 

•  meteorological  conditions 


19 


At  any  particular  geographic  position  the  effect  of  these 
variables  may  or  may  not  be  considered  a  function  of  time,  depending  on 
the  time  periods  of  interest. 

Also  of  importance  is  the  effect  of  these  factors  on  TD  as  a 
function  of  geographic  position.  For  each  different  geographic 
position,  the  electromagnetic  waves  from  each  transmitter  must  traverse 
a  different  route,  with  different  types  of  surface,  topography,  and 
atmospheric  conditions,  to  reach  the  receiver. 

The  effects  of  irregular  terrain  are  practically  time 
invariant,  while  ground  conductivity  and  dielectric  constant  may  vary 
seasonally  (due  to  changes  in  moisture  content  of  the  soil,  changes  in 
vegetation,  or  other  factors) .  Observed  cyclic  annual  variations  on 
the  order  of  0.5  /xsec  are  probably  due  to  seasonal  variations  in  ground 
conductivity  and  dielectric  constant  and  to  seasonally  varying 
atmospheric  conditions. 

Meteorological  conditions  sufficient  to  cause  noticeable 
variations  in  TD  readings  may  change  in  very  short  time  periods  (days 
or  hours) .  Short  term  variations  due  to  meteorological  effects  (cold 
front  and  warm  front  passage  through  the  line  of  propagation)  on  the 
order  of  0.025  jisec  have  been  observed.  These  correlate  with  changes 
in  the  dry  adiabatic  lapse  rate  of  the  index  of  refraction. 
(Ref.l4:p.lll9). 

For  most  differential  Ioran-C  purposes  it  is  necessary  to 
consider  only  those  temporal  propagation  effects  which  occur  in  short 
time  periods  (hours  or  days  at  most,  for  most  applications) .  The 
variation  of  propagation  effects  with  geographic  position  must  also  be 
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considered  for  precise  positioning.  Many  of  these  variations  can  be 
corrected  by  a  computation  of  the  Additional  Secondary  Phase  Factor 
(ASF). 

The  Additional  Secondary  Phase  Factor  (ASF)  is  a  measure  of  the 
systematic  error  resulting  from  differences  in  velocity  of  electromag- 
netic waves  over  ground  with  different  conductivities  and  dielectric 
constants.  The  Defense  Mapping  agency  uses  a  camputerized  system  of 
maps  with  estimates  of  ground  conductivity  (effectively  including 
dielectric  constant)  over  large  areas,  and  calculates  ASF  corrections 
in  the  Loran-C  coverage  areas.  These  corrections  are  published  as 
tables  of  ASF  corrections  to  be  used  in  correcting  loran-C  TD  readings 
from  receivers. 

ASF  corrections  can  be  as  large  as  4  microseconds  (Ref.3). 
The  corrections  provided  in  the  Defense  Mapping  Agency  tables  require 
estimates  of  the  geographic  location  of  the  receiver  to  the  nearest  10' 
of  latitude  and  longitude.  ASF  corrections  based  on  the  known 
positions  of  the  survey  points  used  in  this  experiment  were  provided  by 
DMA  using  a  computer  program  which  calculates  the  corrections  using  the 
Mill ington  Method  (Ref.3:p.V). 

3.  Atmospheric  Noise 

The  effects  of  atmospheric  electrical  noise  are  reduction  of 
the  range  at  which  Loran-C  signals  can  be  effectively  received  and  an 
increase  the  error  in  TD's  cxanputed  by  receivers  in  marginal  reception 
areas.  This  is  primarily  dependent  on  the  circuitry  of  the  receiver. 
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In  marginal  areas  the  effect  of  atmospheric  noise  will  be  an 
increase  in  the  standard  deviation  of  the  TD's  and  increased  difficulty 
in  locking  on  and  maintaining  reception  of  the  weaker  signals. 
Atmospheric  noise  varies  with  time  and  geographic  location,  with  a 
significant  diurnal  variation  (noise  levels  are  higher  at  night) .  In 
this  experiment,  significant  effects  due  to  local  sources  of  electrical 
noise  were  noticed  in  the  Monterey  harbor  area. 

4.  Receiver  Errors 

Receiver  errors  are  of  two  main  types: 

timing  error  caused  by  limitations  of  the  internal  clock 
circuitry  of  the  receiver 
•  errors  caused  by  processing  time  of  the  receiver 

For  a  receiver  at  a  fixed  geographic  location,  processing  time 
is  not  a  source  of  error.  The  receiver  will  be  at  the  same  location 
regardless  of  the  time  required  to  output  the  TD's  or  geographic 
coordinates.  The  main  source  of  receiver  error  at  a  fixed  location 
(other  than  the  noise  considerations  discussed  above)  is  the  timing 
error  in  the  internal  clock  circuitry.  Loran-C  receivers  vary  greatly 
in  their  timing  errors.  The  Racal  Megapulse  Monitor  receiver  used  in 
this  experiment  has  a  specified  timing  error  of  0.02  /xsec  or  better. 

For  mobile  receivers  the  processing  time  of  the  receiver 
circuitry  can  be  important.  After  the  signal  is  received  it  must  be 
processed  and  either  stored  in  the  form  of  electronic  information  or 
visually  displayed.  Error  is  introduced  if  the  receiver  moves  during 
this  time  interval. 
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Obviously,  the  amount  of  error  in  the  determination  of 
geographic  position  from  this  source  depends  on  the  processing  time  and 
the  speed  of  the  receiver.  Receivers  which  output  data  in  the  form  of 
geographic  coordinates  must  calculate  latitude  and  longitude  from  the 
TD  readings  of  the  receiver,  requiring  more  processing  time  and 
resulting  in  greater  error  from  this  source. 

At  the  speeds  normally  used  in  hydrographic  surveying 
(typically  less  than  20  kt)  the  processing  time  is  insignificant  for 
practically  all  receivers.  For  data  which  is  stored  in  the  form  of 
receiver  output  TD's  for  later  correction  and  conversion,  as  with  the 
Racal  Megapulse  receiver  used  here,  error  should  be  insignificant  even 
at  much  higher  speeds. 

D.  Geometric  Dilution  of  Precision 

This  source  of  error  in  geographic  position  results  from  the 
crossing  angle  of  the  loran-C  IDP's  (lines  of  position).  Figure  IV-1 
shows  two  hyperbolic  IDP's  with  standard  deviations  c^  and  o2  (in  time 
difference)  and  crossing  angle  a.  The  optimum  crossing  angle  of  any 
two  IDP's  is  90°,  resulting  in  minimum  error  in  geographic  position  for 
a  given  set  of  standard  deviations. 


23 


E.  Summation  of  Errors 

1.  Absolute  Accuracy  vs.  Precision 

The  absolute  accuracy  of  an  electronic  positioning  system  can 
be  viewed  as  a  measure  of  its  ability  to  give  the  geographic  position 
of  the  receiver.  Precision  is  a  measure  of  the  variation  of  individual 
time  difference  readings.  If  the  receiver  is  set  up  at  a  point  for  a 
sufficient  amount  of  time,  many  of  the  random  errors  can  be  averaged 
out;  the  average  of  a  large  number  of  readings  can  provide  a  basis  for 
determining  absolute  accuracy. 

In  practice  this  is  more  difficult  than  it  would  seem,  since 
some  factors  which  do  not  vary  significantly  over  short  periods  of  time 
(hours  or  days)  may  vary  significantly  over  longer  periods  (weeks  or 
months)  .  It  is  necessary  to  consider  the  time  periods  involved  in 
order  to  determine  which  factors  to  consider  random  and  which  to 
consider  "bias"  in  the  system. 

For  this  experiment,  the  sources  of  random  error  are  considered 
to  be  transmitter  timing  errors,  some  meteorological  conditions, 
atmospheric  noise,  and  most  receiver  errors.  Errors  caused  by  long 
term  "bias"  of  the  system  are  due  primarily  to  ground  conductivity  and 
type  of  surface,  some  meteorological  effects,  and  any  non-random  bias 
in  the  timing  of  the  transmitters. 

Taking  these  factors  into  consideration,  data  taken  at  the 
"known"  point  (Range-7)  of  the  differential  system  was  processed  to 
remove  random  error  by  averaging  over  a  period  of  time,  while  long  term 
error  was  removed  by  comparison  of  the  TD's  to  their  expected  values 
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derived  from  calculations  provided  by  computer  computations  of  the  TD 
at  that  point. 

Errors  which  may  change  within  the  time  frame  of  the  experiment 
were  reduced  by  a  linear  regression  of  the  trend  of  data  taken  at  the 
known  point  and  projected  to  the  correction  at  the  unknown  point 
mathematical ly . 
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IV.   THREE  STATION  AND  DIFFERENTIAL  CORRECTIONS 

A.  Basic  Geometry  of  Loran-C  Hyperbolas 

This  section  provides  a  summary  of  the  information  necessary  for 
the  analysis  of  Loran-C  corrections  for  three  secondary  stations  and 
for  differential  Loran-C  corrections  presented  in  parts  IV. B  and  IV.C, 
adapted  from  Laurilla  (Ref.l)  and  others  as  noted. 

1.   Lanewidth 

The  baseline  of  a  Loran-C  master-secondary  pair  is  the  shortest 
path  over  which  the  transmitted  signals  may  travel  between  the  master 
station  and  the  secondary  station.  On  a  spherical  earth  model  it  is 
the  arc  of  the  great  circle  through  the  positions  of  the  master  and 
secondary  stations.  On  an  ellipsoid  it  is  the  geodesic  from  the 
master  to  the  secondary  station.  In  the  set  of  hyperbolas  of  time 
difference  for  the  master-secondary  pair,  the  separation  between 
hyperbolas  is  a  minimum  constant  value  along  the  baseline. 

A  lane  is  the  distance  between  two  adjacent  hyperbolas  of  unit 
time  difference.  The  unit  normally  used  is  the  i^sec  (10~6  sec) .  In  a 
series  of  hyperbolas,  as  shown  in  Figure  II-l,  it  is  obvious  that  the 
lanewidth  varies  with  location.  At  any  point,  lanewidth  can  be 
viewed  as  the  inverse  of  the  gradient  of  the  ratio  of  the  change  in 
time  difference  to  the  change  in  distance  (Ref.15),  or: 

L  =  {KdTD/dl)!}-1  [IV-1] 


26 


MOST  PROBABLE  POINT 


Figure   IV-1.      Typical   intersection   of   two   lines   of 
position.   Most  probable   point  is   at  the   Intersection  of 

the  two  LOP ' s . 


LOP 


Figure  IV-2.  Typical  intersection  of  three  lines  of 
position.  When  three  or  more  secondary  stations  can  be 
received  at  a  geographic  position,  calculated  position 
depends  on  the  LOP ' s  which  are  used  to  make  the 
calculation. 
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The  lanewidth  at  any  point  in  the  set  of  hyperbolas  can  be 
expressed  in  terms  of  the  lane  width  on  the  baseline.  Baseline  lane 
width  can  be  expressed  by  the  following  equation: 

lb  =  "J"  [IV-2] 

where: 

Lj-,  =  Lanewidth  on  the  baseline  (meters//isec) 

c  =  Speed  of  electromagnetic  waves  on  the  baseline  (m//isec) 
(note  that  time  units  are  in  the  unit  lanewidth) 

Lanewidth  at  any  other  point  (not  on  the  baseline)  can  be 
expressed  by  equation  IV-3a,  where  /?  is  the  angle  subtended  by  the 
baseline  from  the  point.  Note  that  this  angle  is  the  absolute 
difference  between  the  azimuths  from  the  point  to  the  master  and  from 
the  point  to  the  secondary  station. 

L  -  1^(3^03/2)  }-1  [IV-3a] 

or 
L  =  LfcG 
where: 

L  =  Lanewidth  at  a  point  (per  unit  of  time) 
Lj-,  =  Lanewidth  on  the  baseline  (per  unit  of  time) 
/5  =  Angle  subtended  by  the  baseline  from  the  point 
G  =  The  Lane  Expansion  Factor,  G  =  (sin(/3/2)  )-1 

Note  that  the  lane  expansion  factor  depends  only  on  the  angular 
separation  of  the  master  and  the  secondary  stations.   This  obviously 
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depends  on  the  distance  of  the  point  from  the  baseline  center  and  the 
angular  separation  of  the  baseline  and  the  vector  from  the  baseline 
center  to  the  point. 

Equation  IV-3a  provides  a  means  of  calculating  the  normal 
distance  between  any  two  hyperbolas,  provided  they  are  near  enough  that 
their  lanewidths  are  not  significantly  different: 

d  =  Lt  [IV-3b] 

where: 

d  =  distance  in  meters  between  the  two  hyperbolas. 

L  =  lanewidth  at  the  point  of  interest  (in  m/nsec) . 

t  =  time  difference  between  the  two  hyperbolas  in  /Lisec. 

2.  Direction  of  Hyperbola 

The  direction  of  a  hyperbola  at  any  point  on  the  hyperbola  is 
along  the  tangent  at  that  point,  directed  away  from  the  baseline  of  the 
master  and  secondary  stations  used  to  define  the  hyperbola.  For 
applications  in  this  experiment,  Laurilla  (Ref  1)  provides  the  follow- 
ing useful  definition  of  the  direction  of  a  hyperbola:  "...  the 
direction  of  a  hyperbola  at  any  point  coincides  with  the  bisector  of 
the  angle  formed  by  the  lines  joining  the  point  to  the  pair  of 
stations."  (Ref.l:p.92) 

For  purposes  of  this  experiment,  it  is  convenient  to  express 
the  direction  of  the  hyperbola  by  its  azimuth.  This  azimuth  can  be 
expressed  in  terms  of  the  azimuths  (from  north)  of  the  master  and 
secondary  stations  by: 
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Oh  =  (Ofo  +  as)/2    (where  Ici^-as^lSO0)        [IV-4a] 

0^  =  180°  +  (£%  +  oy/2   (where  1 0^-05 1<180°)    [IV-4b] 

where: 

Q^  =  azimuth  of  the  direction  of  the  hyperbola. 

c^  and  as  =  azimuths  of  the  master  and  secondary  from  the  point. 

3.   Applications  of  Geometry  of  Hyperbolas  to  Error  Analysis 

The  geometric  description  of  Loran-C  hyperbolas  in  parts  IV. A.  1 
and  IV. A.  2  is  particularly  important  in  the  analysis  of  error  and  the 
derivation  of  corrections  used  in  this  experiment.  Some  of  the  direct 
applications  are  outlined  in  the  following  paragraphs. 

The  error  expressed  as  an  error  ellipse  reguires  the  directions 
of  hyperbolas  in  order  to  determine  the  orientation  of  the  ellipse. 
The  dimensions  of  the  ellipse  reguire  the  conversion  of  variances  of 
time  differences  of  each  hyperbola  to  variances  of  distance  in  the 
major  and  minor  axis  of  the  ellipse.  Standard  deviations  of  time 
differences  can  be  converted  to  distances  by  the  following  eguations: 

aa  =  atL  [IV-5a] 

ad2  =  at2lj2  [IV-5b] 
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The  correction  for  three  secondary  stations  uses  variances  in 
time  differences  converted  to  distances  as  an  input,   using  equation 
IV-5b.   Converting  time  difference  to  distance  in  this  way  correctly 
weights  the  precision  of  the  IDP's  when  determining  the  correction 
factor  based  on  geometric  considerations. 

B.  Correction  for  Three  Secondary  Stations 

In  situations  where  three  secondary  stations  can  be  received,  the 
result  is  three  LDP's  which  do  not,  in  general,  intersect  at  one  point 
(see  Figure  IV-2)  when  they  are  used  to  calculate  the  geographic 
position,  although  they  must  in  fact  intersect  at  the  common  point  at 
which  the  readings  were  taken  concurrently.  Under  usual  navigation 
conditions  the  two  best  LDP's  are  used  to  determine  a  geographic  fix. 
A  drawback  of  this  method  is  that  it  fails  to  consider  all  of  the 
available  information  (i.e.  one  of  the  available  lines  of  position  is 
not  used  in  the  calculation  of  the  geographic  position) .  By  utilizing 
all  three  LDP's  it  should  be  possible  to  improve  the  accuracy  of  the 
geographic  fix. 

Sections  B.l  and  B.2  will  develop  a  mathematical  adjustment  which 
can  be  used  to  calculate  the  most  probable  corrected  TD  readings  of  a 
Loran-C  receiver  used  as  the  fixed  receiver  in  a  differential  Loran-C 
system  utilizing  three  LDP's  (i.e.  receiving  the  master  and  three 
secondary  stations) .  When  the  correction  is  applied,  the  three 
corrected  TD's  will  intersect  at  the  most  probable  point,  based  on  long 
term  average  values  of  the  TD's.  The  derivation  of  this  adjustment 
could  take  any  of  several  possible  forms,  but  a  primary  consideration 
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here  is  the  need  for  a  form  which  is  easily  applicable  to  computer 
processing.  Section  B.3  describes  a  method  of  applying  the  corrections 
for  non-differential  use,  using  only  one  receiver  and  with  no  known 
fixed  point. 

Section  B.4  will  develop  regression  equations  for  calculating  the 
corrected  TD's  directly  from  time  for  use  as  a  differential  correction 
for  a  remote  receiver.  Note  that  the  methods  developed  in  these 
sections  are  primarily  for  use  in  the  differential  Loran-C  computations 
necessary  for  this  experiment. 

1.  Geometric  Considerations  Used  to  Develop  a  Correction  Factor 

The  following  assumptions  are  made  for  this  derivation: 

•  Each  LDP  has  a  unique  direction  (described  in  Section  IV. A. 2) 

•  Long  term  average  TD's  (designated  TDlav,  TD2av,  and  TD3av)  for 

for  the  known  point  are  available 

•  Variances  of  the  TD  values  (designated  a^2,   a22,  and  a32)  are 

known  or  can  be  calculated 

•  Variations  of  TD's  from  their  average  values  is  in  the  form  of  a 

Gaussian  distribution 

•  Over  the  relatively  small  area  being  considered,  the  surface  can 

be  considered  a  plane  and  LOP's  for  the  same  master-secondary 
pair  can  be  considered  parallel  straight  lines 

Several  things  should  be  noted  about  these  assumptions.  Long 
term  averaging  times  for  this  experiment  will  be  on  the  order  of  a  few 
hours  during  acquisition  of  data;  for  purposes  of  the  derivation, 
standard  deviations  are  converted  between  time  units  (/^sec)  and 
distance  units  (m)  using  Equation  IV-5a;  and,  the  assumptions  of  a 
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FIGURE  IV-3. 


Geometry  used  in  the  derivation  of  a  correction  factor 
for  three  LDP's  not  intersecting  at  the  same  point. 
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Gaussian  distribution,  plane  surface,  and  parallel  LDP's  are  reasonable 
for  this  type  of  data  (Ref  .6:p.8) .  When  two  LDP's  with  different 
standard  deviations  intersect  (see  Figure  IV-1) ,  the  point  of 
intersection  is  the  roost  probable  geographic  position  (Ref .6:p. 8-12) . 
The  probability  of  the  location  of  the  true  position  is  best  described 
using  an  error  ellipse  as  described  in  Section  I1I-B. 

A  typical  situation  is  shown  in  Figure  3,  a  set  of  three 

concurrent  LDP's  intersecting  at  three  different  points.  Point  P  is 
the  location  of  the  monitor  receiver  used  to  determine  differential 
corrections,  with  TDlav,  TD2av,  and  TD3av  being  the  long  term  average 
TD  values  at  that  known  point  (note  that  in  figure  LV-3,  the  LDP's  of 
TDlav,  TD2av,  and  TD3av  are  not  shown,  but  are  parallel  to  TD]_,  TE^, 
and  TD3,  respectively,  and  pass  through  point  P,  as  shown  in  Figure 
TV-4)  .  Although  the  long  term  average  LDP  values  intersect  at  point 
P,  triangle  Pi2pi3p23  represents  a  triangle  of  fixes  formed  by  the  set 
of  TD  values  which  may  occur  from  any  concurrent  set  of  TD's,  which  are 
unlikely  to  equal  the  long  term  average.  Each  set  of  TD  readings 
(which  occurs  in  the  data  every  few  seconds)  will  result  in  a  new 
triangle  of  fixes. 

The  vectors  r^,  r2,  and  r3  represent  the  difference  between  the 
individual  TD  values  in  a  concurrent  set  of  data  and  their  respective 
long  term  TD  values.  The  r  vectors  are  at  right  angles  to  both  the 
long  term  and  current  LDP's.  For  example,  if  the  long  term  average  of 
TD1  (TD1av)  is  29939.200  and  the  value  of  TL\  in  a  set  of  TD  readings 
currently  being  considered  is  29939.300,  the  value  of  r^  for  the 
current  set  of  TD  readings  is  0.100  (all  in  ^tsec)  .   When  converted  to 
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meters  using  equation  IV-3b,  this  is  the  distance  between  the  average 
TD  values  (LOPiav)  and  current  TD  values  (TDjJ  .  Ihese  LDP's  can  be 
considered  as  parallel  straight  lines  over  the  relatively  small 
distances  considered  for  this  correction.  Note  that  there  are  three 
LDP's  in  each  set  of  concurrent  data  representing  TD  values  for  three 
master-secondary  pairs.  The  master  station  is  common  to  each  of  these 
pairs. 

Within  the  triangle  of  fixes  Pi2pi3p23  fQrmed  by  the  three 
LDP's  of  a  set  of  data,  point  P'  represents  the  most  probable  point  of 
the  true  geographic  position  indicated  by  that  particular  concurrent 
set  of  data  compared  to  the  long  term  average  data.  The  distances  of 
P'  from  each  LDP  (the  lengths  r^')  are  related  to  the  relative 
precision  of  each  LDP,  P'  being  closer  to  LDP's  with  higher  precision 
(to  be  discussed  in  more  detail  in  the  next  section) . 

The  problem  must  be  solved  in  general  terms  so  that  any  set  of 
TD's  in  the  data  can  be  corrected  to  a  new  set  of  TD's  which  intersect 
at  the  most  probable  point.  Since  the  triangles  formed  by  all  sets  of 
data  are  similar  triangles,  the  lengths  of  r^',  r2',  and  r3',  will  have 
the  same  proportion  for  each  set.  The  corresponding  r^1  values  between 
sets  of  data  will  be  related  by  a  common  factor,  which  will  be 
designated  as  the  Correction  Factor,  k. 

Other  information  which  should  be  noted  about  Figure  LV-3 
include  the  following: 


The  following  sets  of  lines  can  be  considered  parallel  over  the 
area  represented  by  the  figure:   TDlav  ||  TDX  ;   TD2av  ||  TD2 ; 
TD3av  ||  TD3  ;    rx  ||  rx '  ;    r2  f  r2 '  ;   '  and  r3  f  r3 ' . 


(continued  on  next  page) 
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The  following  sets  of  lines  are  perpendicular:  r^  and  r^'  J_ 
TDlav  and  TC^  ;    r2  and  r2'  J-  TD2av  and  TD2  ;   r3  and  r3'  J_ 
TD3av  and  TD3. 

The  azimuths  (from  north)  of  the  directions  of  the  LDP's 
represented  by  TDlav,  TD2av,  and  TD3av  are  a1;  a2,  and  a3, 
respectively.  These  are  the  same  as  the  azimuths  of  TDlf  TD^ 
and  TD3,  respectively. 

The  following  angles  of  triangle  Pi2pi3p23  are  defined  by: 

a12  =  a2  ~  al  i    a13  =  */2  ~    (a3  ~  al)  ;    aT)d   a23  =  a3  "  al- 
These  angles  will  be  the  same  in  the  similar  triangles  formed 

by  any  lengths  of  r^,  r2,  and  r3. 

LDP's  1  and  3  have  the  largest  absolute  angle  (+  or  -)  between 
them,  and  LDP  2  is  the  intermediate  angle. 


The  next  section  will  provide  a  mathematical  derivation  of  the 
expression  for  a  correction  factor  to  determine  the  lengths  of  r^ ' , 
r2',  and  r3 '  and  the  corresponding  time  adjustments  to  TDlf  1^2*  ar^ 
TE^,  based  on  the  geometry  presented  in  this  section. 

2.  Mathematical  Derivation  of  the  Correction  Factor 

An  important  question  in  this  analysis  is  how  to  weight  the 
relative  precision  of  each  LDP  in  such  a  way  that  it  can  be  used  in 
the  calculation  of  a  correction  factor.  Reference  16  shows  that  for 
this  type  of  data  the  "weights  are  inversely  proportional  to  variance" 
(Ref  .16:p.60) ,  meaning  that  the  lengths  of  the  r^'  corrections  should 
be  proportional  to  the  variances  (in  units  of  length)  of  the 
corresponding  TD  values  for  each  LDP.  Using  units  of  length,  in  the 
general  solution  for  all  sets  of  data,  the  following  equation  applies: 

r^'L^  =  kot2L^2  [TV-6] 

where: 

r^'  =  corrections  for  TD's  1,  2,  and  3  (in  /Ltsec)  . 
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Oj2  =  variances  of  TD's  1,  2,  and  3  (in  /isec2) . 
k   =  the  Correction  Factor  (units  are  meters) 
L^  =  lanewidths  of  TD's  1,  2,  and  3  (in  m//Lisec) . 
Note  that  cjj^L-j2  is  a  unitless  weighting  factor. 

For  any  triangle  there  will  be  some  unique  value  of  k  which 
relates  the  distances  from  point  p'  to  each  LDP  to  the  variances  of  the 
LDP's  according  to  equations  IV-6.  To  obtain  a  general  expression  for 
k  which  can  be  easily  calculated  by  computer,  it  will  be  helpful  to  use 
only  the  known  values  of  the  TD's,  variances,  azimuths,  and  lanewidths 
for  each  LDP.  The  r  values  are  readily  obtained  from  each  set  of  TD 
data  by  subtracting  the  long  term  average  values  from  current  TD 
values.  Values  of  a-j2  and  L^  are  constants  for  this  small  area. 

If  each  LDP  is  considered  a  straight  line  and  a  cartesian 
coordinate  system  is  assumed  with  north  as  the  +y  and  east  as  the  +x 
direction  of  the  coordinate  axes,  then  each  LDP  can  be  specified  by  the 
general  equation  for  a  straight  line: 

y  =  m-jx  +  b^      (i  =  1,  2,  or  3)  [IV-7] 

Three  such  equations  can  be  written  (one  for  each  LDP) ,  giving 
three  linear  equations  in  two  unknowns,  x  and  y.  These  equations  will 
not  have  a  common  solution  when  the  three  LDP's  do  not  intersect  at  a 
common  point. 

For  each  master-secondary  pair,  all  LDP's  of  the  set  will  have 
the  same  slope  (m^)  within  the  small  area  being  considered  (TD^  is 
parallel  to  TDiav  for  any  TD-jJ  .  The  y-intercepts  (the  b^'s)  will  be 
zero  for  the  average  LDP's  intersecting  at  point  P.  Y-intercepts  for 
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other  LDP's  with  non-zero  r  values  can  be  calculated  from  the 
relationships  between  azimuths  (denoted  a^)  and  the  r  values  of  each 
LDP  by  the  following  equation: 


*>i  = 


cos(a^-7r/2)     sin  a^ 


(i  =  1,  2,  or  3)   [IV-8] 
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FIGURE  IV-4.   Typical  average  LDP  through  the  origin  of  a  cartesian 
coordinate  system  and  corresponding  parallel  LDP  of  an  individual 
TD  reading.   Slopes  are  the  same  for  both.   Y  intercept  is  the  b 
value  in  equation  for  the  LDP  (zero  for  the  average  LDP) . 
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Slopes  of  the  TD's  can  be  expressed  in  terms  of  0  as: 

mi  =  -tan  0^  [IV-8b] 

Figure  IV-4  shows  the  assumed  coordinate  system  along  with  the 
LDP's  for  which  the  slopes  and  intercepts  are  described.  When 
corrections  are  applied,  the  corrected  r  values  are  as  follows: 

riCLi  =  riLi  "  toi2Li2  [IV-9] 

Note  that  in  the  equation  for  corrections  (IV-9)  the  units  of 
r^cL^  are  meters.  If  a-j^L-j2  is  considered  a  unitless  weighting  factor, 
then  the  units  of  k  must  be  meters.  This  agrees  with  the  units  of 
Equation  13. 

The  value  of  r^  is  calculated  by: 

ri  =  TDi  -  TDiav  [IV-10] 

where: 

r^  =  deviation  of  an  LDP  from  the  long  term  average  (/xsec) . 

TD^  =  time  difference  reading  of  a  master-secondary  pair  at 
point  P. 

TDiav  =  long  term  average  of  the  TD  for  the  same  pair  at 
point  P. 

Figure  IV-3  shows  that  the  r^'  corrections  must  be  subtracted 
from  r±  and  r3,  and  added  to  r2,  for  the  configuration  shown.  It  can 
be  demonstrated  that  the  r^'  corrections  are  always  of  the  same  sign 
(both  +  or  both  -)  for  the  two  most  widely  spaced  LDP's  and  of  the 
opposite  sign  for  the  intermediate  LOP,  where  the  direction  vectors  at 
point  P  are  in  the  positive  directions  of  increase  of  TD  and  the  angle 
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is  in  absolute  terms  (maximum  being  180°).  This  is  true  for  all 
positive  and  negative  r  values.  Positive  and  negative  k  values  will 
reverse  the  signs  of  all  corrections  in  equations  IV-9  so  that  they 
will  be  correct  for  all  configurations. 

The  corrected  r  values  of  equations  IV-9  can  be  substituted 
into  equations  IV-8,  providing  a  set  of  y-intercepts  (the  b^  values  in 
the  set  of  general  equations  IV-7)  for  the  corrected  IDP's  which 
intersect  at  point  P' .  These  intercepts  contain  the  unknown  Correction 
Factor  k.  Note  that  the  slopes  of  the  corrected  IDP's  are  equal  to 
those  of  the  corresponding  uncorrected  IDP's,  meaning  that  the  m^ 
values  in  equations  IV-7  are  unchanged. 

Substituting  the  corrected  y-intercept  (bjj  values  into 
equations  IV-7  results  in  three  linear  equations  in  three  unknowns:  x, 
y,  and  k: 


rjLi  -  ka^Lj2 

y    =     (-tan  <p±)x  + [IV-lla] 

cos  02 


^1^  +  ^cr22-L22 

y    =     (-tan  02)x  + [IV-llb] 

cos  02 


L3  -  ka32L32 

y    =     (-tan  03)x  +  — [IV-llc] 

cos  03 
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Equations  IV-11  can  be  rewritten  as  the  following  set  of  linear 
equation  in  the  unknowns  x,  y,  and  k: 

al2Ll2  rlLl 

(tan  <p±)x    +  y  +  ( )k  =    — [IV-12a] 

cos  0-^  cos  02. 

a22L22  r2^2 

(tan  02)x  +  Y  ~  ( )k  =     [IV-12b] 

cos  02  cos  02 

a32lj32  r3L3 

(tan  03)x  +  y  +  ( )k  =    [IV-12c] 

cos  03  cos  03 


Solving  equations  IV-12  by  Gaussian  elimination  gives  the  following 
general  expression  for  the  Correction  Factor  k: 


tan  03  -  tan  02       tan  03  -  tan  a^         tan  02  -  tan  0^ 

*ari  _  lrzr2  +  *or3 


cos  0!  cos  02  cos  03 

k  =    

tan  03  -  tan  02        tan  03  -  tan  0^        tan  02  -  tan  0^ 

L^cj-l2 +  l^a22 +  Is?o22 

cos  0^  cos  02  cos  03 

[IV-13] 

where: 

r^  =  differences  from  average  TD's  (/xsec)  . 

k   =  Correction  Factor  for  three  master  secondary  pairs  (m) . 

0i  =  directions  of  increasing  TD  (azimuth  from  north)  . 

Oj2   =  variances  of  LDP's  (part  of  unitless  weighting  factor) . 

L-l  =  lanewidths  (in  meters//Lisec) 
Note  that  L^2^2  terms  in  denominator  are  unitless  weighting  factors. 
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Equation  IV-13  may  be  simplified  by  the  following  substitutions: 

A1  =  L^tan  03  -  tan  02)/ (cos  0X)  [IV-14a] 

A2  =  L2(tan  03  -  tan  01)/(cos  02)  [IV-14b] 

A3  =  L3(tan  02  "  tan  0l)/(cos  #3)  [IV-14c] 

A4  =  Lial2Al  +  L2a22A2  +  L3a32A3  [IV-14d] 

Resulting  in  the  following  equation  for  k: 

k  =   (Airi  -  A2r2  +  A3r3)/A4  [IV-15] 

The  x  and  y  coordinates  of  point  P'  can  be  calculated  by: 


r^cos  02  -  r2Cos  0^  -  ka12cos  0^  -  kCT2  cos  02 

x  =  [IV-16] 

(cos  0]^)  (cos  02)  (tan  02  -  tan  0^) 


y  =  -  (tan  ax)x  +  (rx  -  a^kj/fcos  0X)  [IV-17] 

Two  other  useful  parameters  can  be  calculated  from  the  values 
of  x  and  y.  These  are  the  azimuth  and  range  of  the  of  point  P'  from 
point  P.  These  are  given  by  the  following: 

Azimuth  (from  north) : 

Qppi  =  n/2   -  tan_1(y/x)  (x  >  0)      [IV-18a] 

apPi  =  37T/2  -  tan_1(y/x)  (x  <  0)      [TV-18b] 

Note:  when  x=0  the  azimuth  is  zero  for  positive  y  and  tt   for  negative  y 

Range(P  to  P'):         Rppi  =  (x2  +  y2)^  [TV-19] 
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Using  the  computed  value  of  k  for  a  set  of  concurrent  TD  data, 
a  correction  can  be  made  to  the  TD  values  using  equations  IV-9.  Since 
IV-9  is  in  meters  it  must  be  converted  to  /usee  using  equation  IV-3b. 
The  resulting  corrected  TD's  will  intersect  at  the  point  P' . 

Equations  IV-15,  TV-16,  and  TV-17  do  meet  the  requirements 
stated  at  the  beginning  of  this  section.  The  inputs  are  limited  to 
azimuths,  variances,  long  term  average  TD  values,  and  concurrent  TD 
data,  and  the  values  of  k,  x,  y,  and  the  corrected  TD  values  can  be 
computed  directly  from  these  by  equations  which  can  be  calculated 
easily  by  a  computer. 

3.  Non-Differential  Correction  for  Three  Secondary  Stations 

When  using  the  method  of  Section  IV.  B.  2  with  a  single  Loran-C 
receiver  (i.e.  not  using  differential  methods)  a  correction  can  be 
calculated  based  on  existing  computer  programs.  While  not  used  in  this 
experiment,  the  method  of  correction  will  be  outlined  in  this  section. 

Since  variation  from  a  known  point  cannot  be  used,  the  two  best 
IDP's  can  be  used  as  the  basis  for  this  calculation.  While  not  as  good 
as  using  a  known  fixed  survey  point,  this  will  provide  a  correction 
based  on  three  IDP's.  The  method  is  outlined  below: 


1.  Using  the  two  IDP's  with  the  best  fix  accuracy  (available  from 

existing  computer  programs) ,  calculate  the  geographic 
position  of  the  point  of  intersection.  This  corresponds  to 
point  P. 

2.  From  the  geographic  position  of  P,  calculate  the  TD  of  the 

third  master-secondary  pair  (the  one  not  used  in  the  previous 
step) . 


(continued  on  next  page) 
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3.  Using  point  P  calculate  the  azimuths  (from  north)  of  the  LDP's. 

Number  the  TD's  as  described  in  the  previous  sections,  with 
the  intermediate  TD  as  TD2  and  the  TD's  with  the  greatest 
angular  separation  as  TD^  and  TD3. 

4.  Use  TD  variances  calculated  over  short  time  periods  in  the 

general  area,  or  use  the  values  of  meters/microsecond  cal- 
culated using  equation  TV- 3  as  the  approximate  weights  of 
each  LOP.  Substitute  these  for  the  a2  values  in  equation 
IV-13. 

5.  Assume  r  values  of  zero  for  the  LDP's  with  the  best  fix 

accuracy  used  in  step  one. 

6.  Calculate  the  difference  between  the  observed  third  TD  and  the 

calculated  value  obtained  in  step  2.  Use  this  for  the  third 
r  value. 

7.  Calculate  corrected  r  values  as  in  the  previous  sections. 

8.  Using  any  two  corrected  TD's,  calculate  geographic  position. 

This  is  the  position  corrected  for  three  secondary  stations. 

Note  that  all  of  the  above  calculations  are  based  on  computer 
methods  and  programs  already  in  use  (Ref.  2)  allov/ing  corrected  posit- 
ions to  be  calculated  readily  using  the  method  outlined. 

C.  Differential  Corrections 

One  of  the  main  objectives  of  this  experiment  is  to  apply 
differential  corrections  obtained  at  a  known  point  to  time  difference 
data  taken  at  unknown  points.  One  problem  encountered  was  the 
availability  of  only  one  Loran-C  monitor  receiver,  making  simultaneous 
data  collection  at  two  points  impossible.  Differential  Loran-C  does 
not  require  continuous  updating  of  the  differential  correction  for 
general  navigation  purposes  (hourly  or  daily  corrections  are  possible) , 
but  the  need  for  the  most  accurate  corrections  possible  for  the 
experiment  required  a  continuous  update  of  the  correction. 
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1.  Mathematical  Description  of  the  Differential  Correction 

Data  from  the  experiment  consisted  of  a  series  of  readings  at 
each  point.  Data  was  recorded  at  intervals  of  5.0  seconds  for  much  of 
the  experiment,  although  other  time  intervals  were  also  used.  Data 
consisted  of  time  versus  TD  for  up  to  four  possible  master-secondary 
pairs  and  the  signal  to  noise  ratios  for  each,  in  the  following  format: 

Time  (GMT)    GRI   TT^       TD2   TD3   TD4   S/U±       S/N2   S/N3   S/N4 

Only  three  master-secondary  pairs  were  receivable  in  the  area 
of  the  experiment.  One  survey  point  (Range  7}  was  designated  the 
"known"  point,  from  which  differential  corrections  were  calculated  for 
data  taken  at  the  other  points.  Over  a  period  of  several  days,  data 
was  collected  at  the  known  point  for  periods  of  20  to  30  minutes  at 
irregular  intervals  spaced  at  about  one  to  three  hours  (data  was 
collected  at  the  other  points  in  the  intervening  periods  as  described 
in  Section  A. 2) . 

To  provide  estimates  of  the  continuous  TD  readings  at  the  fixed 
point  during  all  time  periods  (including  the  periods  between  data 
collection  at  that  point) ,  a  least  squares  curve  fitting  process  is 
used.   The  result  is  three  equations:   TD  as  a  function  of  time,  for 
each  master-secondary  pair. 

These  three  equations  can  be  used  to  provide  an  estimate  of  the 
TD  values  at  the  fixed  point  at  any  time.  The  difference  between  the 
calculated  TD  at  a  certain  time  and  the  long  term  average  TD  is  the 
Differential  Correction  at  that  time.  The  ASF  corrections  at  the  fixed 
point  and  an  estimate  of  the  ASF  correction  at  the  unknown  point  (as 
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described  in  Section  II-B)  also  enter  into  the  Differential  Correction 
at  the  unknown  point.  An  equation  is  written  for  each  LDP: 

Cd  =  TDC  -  TDav  +  (ASFU  -  ASFf)  [IV-20] 

where: 

C$    =  Differential  Correction  at  the  unknown  point. 

TDC  =  TD  (at  known  point)  calculated  from  regression  equation. 

TD^y  =  Long  term  average  TD  at  the  known  point. 

ASF  =  ASF  correction  at  the  known  and  unknown  points. 

2.  Least  Squares  Regression  Equations  for  Time  Difference 

Linear  regression  equations  for  each  TD  can  be  written  in  the 
following  form  (Ref.17): 

TD  =  b0  +  b]t  [IV-21] 

The  b^  coefficients  of  the  regression  equation  can  be  found  by 
solving  the  following  set  of  linear  equations: 


nb0     +  ^Et^   =  STD-l 

[IV-22a] 
boEtj^   +  b^Et^2  =  Et-jTDi 


An  alternate  method  of  calculating  the  b^  values  directly  is: 


bx  = 

b0  = 

nEtjTD-L  -  (EtjJ  (ETDj^) 

nEti2  -  (Eti)2 
ETD-l       ETD^ 

V-i 

-  b± 

[LV-22b] 


n         n 
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In  equations  IV-22,  Et-^,  St^TD^,  etc.,  refers  to  the  sum  of 
the  time  values,  the  sum  of  the  products  of  time  and  time  difference 
values,  etc. ,  for  the  whole  set  of  data  available  at  the  fixed  point. 
The  total  number  of  sets  of  TD  vs  t  data  is  n. 
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V.      RESULTS  AND  OCNCLUSICNS 


Part  A  of  this  chapter  presents  the  results  of  calculations  based  on 
the  theoretical  principles  previously  discussed  and  data  obtained  from  the 
experiment.  Part  B  provides  a  summary  of  the  results  obtained,  how  these 
results  differ  from  what  would  be  expected  based  on  the  theory,  possible 
reasons  for  any  variations  from  expected  outcome,  and  a  general  conclusion 

A.   Calculations 

1.  Azimuths,  Directions  of  Hyperbolas,  and  Directions  of  TD's 

Azimuths  of  the  Loran-C  stations  from  the  survey  points  are  one 
of  the  basic  requirements  of  several  calculations  used  in  this 
experiment.  Station  positions  of  Range-7  and  Luces  Point  (the  primary 
stations  used  in  this  experiment)  and  other  points  were  obtained  from 
Horizontal  Control  Data  published  by  the  U.S.  Dept.  of  Commerce,  Coast 
and  Geodetic  Survey  (see  Appendix  A) .  Station  positions  of  the  Loran-C 
transmitters  were  obtained  from  Loran-C  Data  Sheets  provided  by  the 
Defense  Mapping  Agency,  Hydrographic/Topographic  Center  (see  Appendix 
B) .  For  convenience  in  avoiding  datum  conversion,  the  North  American 
Datum  1927  (NAD  27)  was  used  because  published  positions  based  on  this 
datum  were  available  for  all  of  the  station  positions  used  in  the 
experiment.  Station  positions  are  summarized  in  Table  V-l. 
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Azimuths:  Azimuths  from  Range-7  and  Luces  Point  to  Loran-C 
stations  9940-M,  W,  X,  and  Y  were  calculated  using  a  Fortran  77 
subroutine  provided  by  DMAHTC  (Appendix  C).  Operation  of  the, 
program  was  checked  using  known  data  points  and  azimuths.  Results  of 
these  computations  are  presented  in  tabular  form  in  Table  V-2.  For 
convenience  in  visualizing  the  angular  relation-  ships,  results  are 
also  presented  in  graphic  form  in  Figures  V-la  and  V-2a. 

Directions  of  Hyperbolas:  The  directions  of  hyperbolas  for 
9940-W,  X,  and  Y  were  calculated  using  equations  IV-4  and  the  azimuths 
obtained  as  described  in  the  previous  paragraph.  A  sample  calculation 
follows: 

Hyperbola  Direction  for  9940-X  at  Range-7: 

(Azimuth  of  9940-X)  +  (Azimuth  of  9940-M) 
HDr7  =  

2 

343.071686  +  38.358551 

=  190.715118° 

2 

(from  North) 

Decimal  degrees  are  used  in  the  calculations  for  convenience. 
Results  of  the  calculations  are  again  presented  in  Figures  V-lb  and  2b. 

Direction  of  Increase  of  TD:  Directions  of  increasing  TD  are 
at  right  angles  to  the  hyperbola  directions  and  increase  in  the 
direction  of  hyperbolas  which  approach  the  master  station.  This 
involves  either  adding  or  subtracting  90°  from  the  direction  of  the 
hyperbolas.  For  checking  purposes,  a  Loran-C  chart  of  the  area  was 
used.  Results  are  again  presented  in  Figures  V-l  and  2. 
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TABLE  V-l.   Station  Positions  (NAD  27)  . 


Station 

Latitude 

Longitude 

Range  7 

36° 

391 

02.47787"  N 

121° 

49'  08.58202"  W 

Luces  Point 

36° 

38' 

10.524"  N 

121° 

55'  38.399"  W 

9940-M 

39° 

33' 

07.046"  N 

118° 

49'  52.241"  W 

994 0-W 

47° 

03' 

48.594"  N 

119° 

49'  52.241"  W 

9940-X 

38° 

46' 

57.472"  N 

122° 

29'  40.050"  W 

9940-Y 

35° 

19' 

18.342"  N 

114° 

48 '  13.946"  W 

TABLE  V-2.  Aziinuths  From  Range-7  and  Luces  Point. 


Range-7  to 

Degrees , min , sec . 

Decimal  Degrees 

9940-M 
994 0-W 
9940-X 
9940-Y 

38°  211  30.784" 

7°  32'  22.823" 

343°  04'  18.070" 

101°  47'  36.229" 

38.358551° 

7.539673° 

343.071686° 

101.793397° 

Luces  Point  to 

38°  51'  08.262" 

7°  45'  43.837" 

344°  23'  49.981" 

101°  361  16.319" 

38.852295° 

7.762177° 

344.397217° 

101.604533° 

9940-M 
994 0-W 
9940-X 
9940-Y 
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9940-W  7.539673* 


9940-X  343.071686° 


9940-M  38.358551* 


9940-Y  101.793397 


FIGURE  V-la.  Azimuths  Fran  Range-7  to  9940-M,W.X,and  Y. 


9940-Y  250.075974 


9940-W  202.949112' 


// 


9940-X     190.715118' 


FIGURE  V-lb.      Directions  of  Hyperbolas  at  Range-7. 


9940-Y     340.075974* 


9940-X     100.715118* 

9940-W     112.949112* 
FIGURE  V-lc.     Directions  of  Maximum  Increase   (Gradient)   of  ID'S  at  R-7. 
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9940-W     7.762177' 


9940-X     344.397217° 


9940-M     38.852295* 


9940-Y     101.604533' 


FIGURE  V-2a.     Azimuths  Fran  luces  Point  to  9940-M,W.X,and  Y. 


9940-Y     250.228424" 


9940-W     203.307236°  / 


9940-X      191.624756° 


FIGURE  V-2b.     Directions  of  Hypertolas  at  Iuoes  Point. 


9940-Y      340.228424 


9940-X      101.624756° 


9940-W     113.307236° 
FIGURE  V-2c.      Directions  of  Maximum  Increase   (Gradient)   of  TO's  at  IP. 
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2.  Lanewidths 

Lanewidths  were  calculated  from  Equation  IV-3a  using  the 
azimuths  computed  in  Section  V-A-l.  Note  that  the  units  of  lanewidth 
are  m/nsec.  With  a  standard  lane  of  1  iisec  the  width  of  a  lane  is  the 
same  number  in  meters.  Velocity  of  transmission  was  calculated  using  c 
and  the  index  of  refraction  specified  in  the  table  of  Loran-C  constants 

supplied  by  the  Defense  Mapping  Agency  (Ref.  4).  Sample  calculations  of 
lanewidth  follow: 


c 

Velocity  of  Transmission  =  [V-l] 

n 


2.99792453  x  108  m/sec 


1.000338 


=  2.996911624  x  108  m/sec 


Lanewidth  of  994 0-Y  at  Range-7  for  standard  1  /xsec  lane: 
(10~6sec) (2.996911624  x  108  m/sec) 


^7 


101.793397°  -  38.358551° 
2  sin  


=  285.024  m 


Results  of  the  calculation  of  lanewidths  are  given  as  the  L 
values  in  Table  V-3.  In  the  table,  L^,  1%,  and  L3  correspond  to  the 
lanewidths  for  9940-W,  X,  and  Y,  respectively. 
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3.  Average  TD's,  Standard  Deviations,  Variances 

Because  of  the  large  amount  of  data  involved,  it  was  necessary 
to  compute  statistical  information  using  a  computer  program.  Since 
data  was  in  the  form  of  three  TD's  (one  for  each  master-secondary  pair) 
every  five  seconds,  and  readings  were  taken  for  periods  of  20  to  40 
minutes  at  each  station,  a  Fortran  77  program  with  subroutines  for 
statistical  information  was  written  (see  Appendix  C) . 

Station  Range-7  was  the  "known"  station  of  the  differential 
system,  at  which  more  readings  were  taken  than  any  other.  For  this 
reason,  and  because  regression  equations  of  TD  vs  Time  were  needed  only 
at  that  station,  statistical  information  was  calculated  only  for 
station  Range-7. 

Mean  of  the  TD's  was  computed  by  the  following  (Ref.  18): 


ETDi 

Mean  TD  -  [V-2] 

n 


Variances  of  the  TD's  were  calculated  using  the  following  (Ref.  18) 


!•  2  _    fvi'n-  ^  2 , 


ZTDj^   -    (ZTDjJVN 
N  -   1 


s2  =  —  [V-3] 


Standard  deviations  were  calculated  from  variances  by  (Ref.  18): 
formula  (Ref. 18): 

s  =   (s2)^  tv-4] 
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Equations  V-2,  V-3,  and  V-4  were  incorporated  into  a  Fortran  77 
computer  program  which  was  used  to  calculate  statistics  of  TD's  for 
9940-W,  X,  and  Y. 

It  should  be  noted  that  these  statistics  were  in  terms  of  time 
differences,  while  in  some  cases  it  was  necessary  to  know  standard 
deviations  and  variances  in  terms  of  distances.  Equations  IV-5  were 
used  for  this  conversion  where  necessary. 

Corrected  vs  Uncorrected  Data;  One  use  of  the  statistics  of 
the  data  was  to  compare  the  standard  deviation  of  raw  data  to  the 
standard  deviation  of  data  which  was  corrected  for  three  secondary 
stations.  It  was  considered  unlikely  that  the  three  station  correction 
could  remove  bias  from  the  TD  readings  (error  in  terms  of  absolute 
position)  because  most  of  this  is  due  to  systematic  error  rather  than 
random  variations.  It  was  believed  that  the  three  station  correction 
would  be  most  effective  in  reducing  the  effect  of  random  variations  in 
time  differences.  The  statistics  on  uncorrected  and  corrected  data  show 
results  consistent  with  this  expectation  (Table  V-7). 

4.  Error  Ellipse  and  Circular  Error  at  Range-7 

Range-7  was  used  as  the  known  point  of  the  differential  system, 
and  a  larger  amount  of  data  was  collected  at  that  point.  Using  the 
crossing  angles,  standard  deviations  and  variances  of  9940-X  and  Y  (the 
best  pair  at  that  point) ,  dimensions  and  orientation  of  the  standard 
error  ellipse  and  CEP  50%  and  90%  circle  radii  were  calculated  for  the 
point  using  the  information  outlined  in  section  III-B. 
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The  standard  error  ellipse  and  error  circles  for  50%  and  90% 
were  plotted  over  a  scatter  plot  of  data  from  Range-7.  The  result  is 
shown  in  Figure  V-3. 

Calculations  of  the  standard  error  ellipse  dimensions  were  as 
follows : 

Using  9940-X,Y: 

°1     =     a9940-Y  =   (0.072546) (285.024)   =  20.677m 

°2     =     a9940-X  =   (0.068638) (322.024)   =  22.103m 

<p     =  250.075974°  -  190.715118°   =  59.360856° 

Calculating  the  semi-major  and  semi-minor  axes  of  the  ellipse: 

20 . 6772+22 . 1032+ ( (20 . 6772+22 . 1032) 2-4 (20 . 677222 . lO32sin20) ) * 


°x2     = 


ox       =  30.595  m 


2sin20 


20 . 6772+22 . 1032- ( (20. 6772+22 . 1032) 2-4 (20. 677222 . lO32sin20) ) * 

°y2    =   


2sin20 


Oy       =     17.362  m 


Orientation  of  the  major  axis  of  the  ellipse: 

20.6772sin  2(59.360856°) 

tan  29  =  

20.6772  cos  2(59.360856°)  +  22.1032 

9  =  26.474268° 

Azimuth  of  major  axis  =  250.075974°  -  26.474268  =  223.601716° 

=  223.601716°  -  180°  =  43.601716°  (from  North) 
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N 


100 


Figure  V-3.   Scatter  Diagram,  Standard  Error  Ellipse,  and  CEP50  and 
CEP90  Error  Circles,  for  Range-7  Data,   (see  text  for  details) 
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Using  tables  in  Reference  6,  pp.  18-20,  radii  of  the  50%  and 
90%  error  circles  were  calculated  as: 

CEP50  =  27.936  m  CEP90  =  63.226  m 

Because  the  data  included  in  the  plot  of  Figure  V-3  includes 
two  intervals  of  data  collection,  the  points  are  not  as  evenly 
distributed  as  they  would  otherwise  be.  For  this  reason  there  are 
fewer  points  near  the  center  of  the  plot  than  would  normally  be  the 
case.  The  error  circles  and  ellipse  seem  to  be  reasonably  correct  for 
the  given  data. 

5.  Correction  Factor  for  Three  Secondary  Stations 

The  correction  for  three  secondary  stations  was  designed  as  a 
technique  to  reduce  the  random  error  associated  with  time  difference 
data.  It  was  expected  to  increase  the  precision  of  the  position  rather 
than  the  absolute  position. 

For  a  variety  of  reasons  there  will  be  some  covariance  between 
concurrent  time  difference  readings  from  different  master-secondary 
pairs.  The  correction  for  three  secondary  stations  is  not  expected  to 
reduce  the  error  associated  with  this  type  of  random  error.  The  non- 
correlated  random  error  in  TD  readings  should  be  reduced  by  this 
correction  to  some  extent. 

In  order  to  reduce  the  number  of  computations  and  the  time 
required  to  process  data,  most  constant  factors  associated  with  the 
correction  factor  can  be  computed  in  advance  and  used  as  constants  in 
the  computer  program.   This  reduces  the  correction  to  only  a  few  lines 
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of  relatively  simple  conputation.  Because  of  the  large  amount  of  data, 
this  reduces  significantly  the  time  required  to  compute  corrections. 

In  the  calculations  of  adjustments  to  positions  of  Range-7  and 
Luces  Point,  the  only  factors  which  differ  significantly  are  the  long 
term  time  differences  at  the  two  points.  In  this  case  the  same  A^ 
factors  can  be  used  for  correcting  both  positions  for  three  secondary 
stations.  In  cases  where  some  of  the  parameters  may  differ 
significantly,  it  would  be  necessary  to  calculate  different  A^  factors 
for  each  position.  Table  V-3  gives  parameters  for  input  to  the  Fortran 
77  correction  program  for  calculation  of  the  A±  parameters  and 
calculation  of  the  corrections  to  data. 


TABLE  V-3.    Parameters  used  in  calculation  of  correction 

factor  k . 


Parameter 

Value 

01 

112.949112° 

02 

100.715118° 

03 

340.075974° 

°1 

0.150606  usee 

°2 

0.068638  /xsec 

°3 

0.072546  fisec 

Ll 

563.934  m/Msec 

^2 

322.966  m/Msec 

L3 

285.024  m/fiSGC 
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6.  Regression  Equations  for  Time  Difference  vs  Time 

Since  only  one  Loran-C  monitor  receiver  was  available,  it  was 
necessary  to  develop  a  method  of  estimating  TD's  at  the  known  point 
during  times  when  differential  corrections  were  required  at  the  unknown 
points.  Linear  regression  equations  of  TD  vs  time  were  used  for  this 
purpose. 

Using  equations  IV-21  and  TV-22,  the  data  from  Range-7  was 
processed  to  obtain  the  required  summations  of  time  and  time  difference 
necessary  to  calculate  the  bo  and  b^  coefficients  of  the  linear 
regression  equation.  Only  Range-7  data  from  the  nearest  time  periods 
before  and  after  the  differential  correction  period  for  Luces  Point 
data  were  used. 

The  procedure  for  processing  the  data  was  first  to  apply  the 
correction  for  three  secondary  stations  to  the  Range-7  and  Luces  Point 
data  during  the  time  periods  of  interest.  Range-7  data  was  then 
processed  to  obtain  the  necessary  summations  for  the  calculation  of  the 
b  parameters  of  the  regression  equations.  Using  the  regression 
equations  for  Range-7  data  and  long  term  TD  averages  at  the  two  points, 
differential  corrections  for  the  Luces  Point  data  were  calculated  as 
described  (Section  V-7).  Fortran  subroutines  (Appendix  C)  were  added 
added  to  the  statistics  program  to  obtain  the  summations  and  a  separate 
program  was  used  for  the  differential  correction  of  Luces  Point  data. 
Results  of  the  corrections  were  compared  to  calculated  values  of  TD's 
provided  by  the  Defense  Mapping  Agency,  Hydrographic/Topographic 
Center.  TD  vs  time  is  plotted  in  Figures  V-4  along  with  the  regression 
line  for  the  data.  Regression  equations  are  in  Table  V-4. 
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Figure  V-4a.   TD  vs  time  and  the  regression  line  for  9940-W 
for  Range-7   data  corrected  for  three  secondary  stations. 
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Figure  V-4b.   TD  vs  time  and  the  regression  line  for  9940-X 
for  Range-7  data  corrected  for  three  secondary  stations. 
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Figure  V-4c.   TD  vs  time  and  the  regression  line  for  9940-Y 
for  Range-7  data  corrected  for  three  secondary  stations. 
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TABLE  V-4.  Regression  Equations  of  TD  vs  Time  at  Range-7, 


994 0-W: 

T°W 

= 

16316.51564   -  0.00807243t 

994 0-X: 

TDx 

= 

27522.33725  -  0.01385199t 

994 0-Y: 

TDy 

= 

42749.29371  +  0.01056707t 

time  units:  hours  with  minutes  and  seconds  in  decimal  hours. 
7.  Differential  Corrections 

Equation  IV-20  describes  the  differential  correction  which  is 
applied  to  TD  readings  at  an  unknown  point  based  on  TD's  at  a  known 
point.  In  this  experiment  the  positions  of  all  points  are  known,  based 
on  their  published  positions  from  survey  data  (see  appendices)  .  One 
point,  Range-7,  was  chosen  as  the  known  point  for  purposes  of 
differential  corrections.  Data  from  one  other  point  was  chosen  as  the 
unknown  point  for  purposes  of  comparison  to  its  known  position. 

Three  things  must  be  known  to  calculate  a  differential 
correction  based  on  equation  IV-20.  Rather  than  use  long  term  averages 
in  this  equation,  the  expected  TD's  at  both  points  were  computed  by  a 
standard  Loran-C  program.  Current  TD's  at  Range-7  were  calculated  from 
the  regression  equations  described  in  the  last  section.  The  ASF 
corrections  in  equation  IV-20  were  calculated  for  each  point  by  the 
Defense  Mapping  Agency  using  a  computer  program  which  is  used  in  the 
preparation  of  ASF  correction  tables,  but  which  can  be  used  for 
determining  the  correction  at  individual  points.   The  difference  in  ASF 
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corrections  is  used  in  equation  IV-20.   Table  V-5  gives  expected  TD 
values  (used  as  long  term  averages)  at  the  two  points  of  interest: 

TABLE  V-5.   Calculated  TD's  used  as  long  term  averages  (in  /xsec) . 


Range-7 

Luces  Point 

9940-W 

16314.68 

16300.09 

9940-X 

27523.13 

27496.79 

994 0-Y 

42749.49 

42755.87 

ASF  corrections  were  approximately  the  same  at  both  points  due 
to  the  nearness  of  the  points.  Differences  in  ASF  corrections  used  in 
equation  IV-20  were  very  small:  0.005,  0.002,  and  0.003  /isec  for  9940- 
W,  X,  and  Y,  respectively.  Corrections  of  this  magnitude  are  probably 
less  than  the  accuracy  of  the  data  used  to  calculate  them,  but  at 
longer  ranges  these  may  be  large  enough  to  be  significant. 

The  final  equation  for  the  differential  corrections  based  on 
TD's  at  Pange-7  are  in  Table  V-6.  Note  that  t  is  Pacific  Standard  Time 
in  hours,  with  minutes  and  seconds  expressed  as  a  decimal. 


TABLE  V-6.  Differential  Corrections  (C^  in  /xsec,  t  in  hours  PST) 


^dw 

= 

-1.83064  +  0.00807243t 

°dx 

= 

0.79575  +  0.01385199t 

Cdy 

= 

0.19829   -   0.01056707t 
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8.  Accuracy  of  Position  Fix 

Figure  V-5a  shows  a  plot  of  uncorrected  Range-7  data  based  on 
expected  TD's  calculated  from  the  known  geographic  position.  The  grid 
is  in  meters,  oriented  north-south  and  east-west,  and  with  the  center 
at  the  published  geographic  position  of  Range-7.  When  the  Correction 
For  Three  Secondary  Stations  is  applied  and  the  average  of  the  Range-7 
data  is  placed  at  the  origin,  the  result  is  as  shown  in  Figure  V-5b. 

Note  that  the  plot  is  now  centered  on  the  origin  due  to  the 
fact  that  all  absolute  error  (bias)  has  been  removed  from  the  data. 
The  three  station  correction  has  also  decreased  the  standard  deviation 
and  variance  of  the  Range-7  data  as  shown  in  Table  V-7. 

Figure  V-6a  shows  luces  Point  data  before  application  of  any 
corrections.  These  points  are  plotted  based  on  the  expected  TD's  at 
the  known  coordinates  of  the  point.  The  center  of  the  uncorrected 
Duces  Point  data  is  at  x  =  375.820  m,  y  =  85.501  m,  and  at  a  range  of 
385.423  m  and  at  an  azimuth  of  77.183°  (from  North)  from  the  origin. 
Radii  of  the  CEP50  and  CEP90  error  circles  are  14.9  m  and  28.4  m, 
respect  ively . 

Figure  V-6b  shows  Luces  Point  data  after  the  application  of 
both  three  station  and  differential  corrections.  Again,  the  origin  is 
at  the  known  position  of  the  point.  For  visual  comparison  purposes, 
the  same  scale  is  used  for  both  Figures  V-6a  and  b.  For  the  corrected 
data,  the  center  of  the  data  is  at  x  =  47.213  m  and  y  =  -7.015  m,  with 
a  range  of  47.731  m  and  at  an  azimuth  of  98.451°  (from  North)  from  the 
origin.   C£P50  and  CEP90  radii  were  12.6  m  and  26.0  m,  respectively. 
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Figure  V-5a. 


Scatter  Diagram  For  Range-7  Uncorrected  Data.   Scale  is 
in  meters  oriented  N-S  and  E-W,  centered  on  the  time 
differences  calculated  from  the  published  coordinates 
of  the  point. 
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Figure  V-5b. 


Scatter  Diagram  for  Range-7  Data  Corrected  for  Three 
Secondary  Stations  and  Differential.   Scale  is  in 
meters  oriented  N-S  and  E-W.   Center  of  data  is 
corrected  to  correspond  with  calculated  time 
differences  for  the  published  coordinates  of  the  point. 
See  text  for  details. 
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Figure  V-6a.   Scatter  Diagram  For  Luces  Point  Uncorrected  Data.  Scale 
is  meters  oriented  N-S  and  E-W,  centered  on  the  time 
differences  calculated  from  the  published  coordinates 
of  the  point. 
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Figure  V-6fc>. 


Scatter  Diagram  far  Dices  Point  Data  With  Differential 
Corrections  Based  on  Range-7  Data  and  With  Correction 
For  Three  Secondary  Stations.  Origin  of  plot  is  at  the 
time  differences  calculated  from  the  published 
coordinates  of  the  position. 
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TABLE  V-7.   Standard  Deviation  and  Variance  of  Range-7  Data  With  and 
Without  the  Correction  Poor  Three  Secondary  Stations. 


Pair 

Without  Correction 

With  Correction 

a 
9940-X 

0.068638 
0.004711 

0.051161 
0.002617 

a 
9940-Y 

o2 

0.072546 
0.005263 

0.052081 
0.002712 

B.  Conclusions 

Based  on  the  results  given  in  part  A,  the  absolute  accuracy  of 
the  Luces  Point  Data  was  improved  from  an  average  error  of  385.423 
meters  to  an  average  error  of  47.731  meters.  Practically  all  of  this 
improvement  was  a  result  of  application  of  the  differential  correction. 

Precision  of  the  Luces  Point  data  was  improved  from  14.9  m  and 
28.4  m  for  the  50%  and  90%  CEP  radii,  uncorrected,  to  12.6  m  and  26.0  m 
for  the  corresponding  radii,  corrected.  Most  of  this  improvement  is  a 
result  of  the  correction  for  three  secondary  stations. 

These  figures  are  based  on  the  processing  of  486  sets  of  TD 
data  from  Range-7,  the  "known"  point  of  the  differential  system,  and 
252  sets  of  TD  data  at  Luces  Point.  Each  set  of  data  included  time  and 
three  TD's,  one  for  each  master-secondary  pair  (9940-W,  X,  and  Y) . 

Under  ideal  conditions  a  better  absolute  correction  would  be 
expected  for  these  points.  From  an  examination  of  the  characteristics 
of  the  data,  the  weakest  point  is  probably  the  lack  of  a  continuous  set 
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of  TD  readings  at  Range-7.  Error  in  the  estimate  of  TD's  based  on  the 
linear  regression  lines  may  be  as  much,  or  more  than,  the  48 -meter 
average  error  in  the  corrected  positions  calculated  for  luces  Point. 

Other  sources  of  error  may  include  the  fact  that  some  of  the 
loran-C  signals  traversed  some  seawater  path  near  Luces  Point,  with 
unknown  results  on  TD  readings  at  that  point.  This  type  of  error  is 
known  to  result  under  similar  conditions,  but  no  means  of  estimating 
its  magnitude  is  available  at  this  time. 

In  summary,  the  absolute  positions  calculated  from  differential 
corrections  of  Loran-C  data  seem  to  be  a  significant  improvement  over 
the  positions  calculated  from  raw  data.  While  the  correction  for  three 
secondary  stations  does  not  add  much  to  the  estimate  of  positions  in 
this  experiment,  the  improvement  of  50%  CEP  radii  by  about  15%  does 
seem  significant  enough  to  warrant  further  consideration. 
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C        THIS  PROC-RAH  READS  0*1*  FROn  A  FILE  OF  TIFlE.  A.  I.     C 

C        CONTINUES  TO  THE  END  OF  THE  FILE,  CALCULATING  THE  HEAN, 

C        VARIANCE.  AND  STANDARD  DEVIATION  OF  THE  A,  B,  AND  C  COLUMNS. 

C        AND  PRINTS  OUT  THE  RESULTS  INTO  A  FILE  tdout.f 

double  precision  t lie . t iae 1 , t iee2 , t iae3, e, b, c . na, ne. nb, nc , 

«>Ut,  lull.  SUaC  .  inn.  lo«  a  2  ,  Suab2 ,  Suac2.  I  .»0  ',  ,  lult]  ,  (UtCl   6  1  ,  b]  ,  C  1 

open  (unit'l.file  «  'out9.f',  for*  *  'foraatted'. 
♦access='seouentlel'  ) 
C       open  (unit«2,file  »  'out6.f',  fora  «  'foraatted', 
C      «eccess« ' sequent ial '  ) 

rewind  1 
C       reedl 1 , '  )  frst In 
C       renal  1 , *  )  etc  In 

SO  write  (*,*)  t iael , tiae2, t iaeJ.al ,bl ,cl,na, na.nb, nc 

tiae  •  tlael  ♦  tiae2/60.  «  tiee3/3600. 
C       write  (  2, I  I  t  1  ae , a  1 ,  b  1 , c  1 

sua  ■  sua  ♦  1 . 


aura  '    suae  « 

a 

suaei 

«  suasl 

♦ 

al 

iuaa2 

»  lur«; 

4 

e"e 

suab  i 

i  suab  < 

b 

suabl 

»  suat : 

♦ 

bl 

SUBb2 

=  suab2 

♦ 

b'b 

suae  i 

■  suae  ♦ 

c 

tunc  1 

=  sumel 

♦ 

cl 

suae  2 

'     suae  2 

♦ 

c  "c 

reed(l,2,err=10,iostat»ioerr)  tiael,ti*e2,tiaeJ,el.a,bl,b. 

♦  el,c 

it     (ioerr)    10. SO, 10 
10    suit    =     suit    -     1  . 
wr  1 1  e  ( • . 6 )    ioerr 
call    s2(su»e2.suaa,sui. ,  s  u  a  a  1  ) 
cell    s2fsu»b2, suab, sua, suabl) 
cell    s2  I  sui»c2.  suae  ,  sua,  suae  1  ) 
2    f oraet (t 1 . d2. 0, tt , d2.0, t 7, d2.0, tie, d8. 3, t21,dS.3, t27, 

♦08. 3, t JO. dt. 3, t36.de    3, t39,dS. 3) 
i    foraetlf3.0.2i,f3.0,2«.f3.0,2i,fl0.3,2i,fl0.3,2i,fl0.3.2i, 

♦  U,  it,  ia.laj 

t    f oraet ( ' i  oerr    *    ' , i  J ) 

8  for»et(fl2.8.2i,fl0.3.2i.fl0.3,2i,fl0.J) 

9  for  mat (' dot e    error') 
40    stop 

end 

subroutine  s2 ( ■ 2 , i , d , s 1 ) 

double  precision  i  2 , i , d, aean, var ian . sd, si 

aeer,  =  sl/d 

verier.  *  (  i  2- i  "  « /d  )  /  I  d- 1  ) 

sd  *  dsor t ( ver  ian ) 

write!", 12)  aeen,  ver i an, sd 

wrlte(2,12)  aean, ver 1  en, sd 
ret  urn 
12  foraatl'aean  ■  ', fl 2. 3, 2> ,' variance  *  ' , f 1 2. 8. 2i , ' sd  =  ',fl2.8) 
end 
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10    re.  

20    re.  --- 

50    re.  

40    rc,  

50    re.    

60  re. 
70  open  "1" 
BO  Input  II 
90  r  »  row* 
95  re. 
100  re.  --- 
105  re. 
110  for  j  » 
120  Tor  k  • 
150  input  I 
140  if  ►  () 
150  Alfj.k) 
160  ne>t  k 
170  neit  j 
175  re. 
180  re.  --- 


PROGRAM  FOR  tAUSSIAN  ELIMINATION  WITH  riVOTAL  CONDENSATION 

FOR  SOLUTION  OF  SYSTEMS  OF  LINEAR  EOUATIONS   

- WITH  UP  TO  10  UNKNOUNS    


II,  "rtest.f  :  dl.  All  10.  1  1  )  :du  BIMll.di.  Hl(ll):di.  Hill) 

r ovs ! i  coll! 

:  c  »  coll!  :  rowl  «  1  :  coll  «  1  :  highrow  •  1 


Input  the  ■■tri i : 


1  to  rows ! 
1  to  coll! 
1,  All  j.k) 
cols!  then  160 
»  -1  "  Alfj.k) 


begin  nin  progre. 


185  re. 

190  gosub  500 

200  gosub  890 

210  eosub  650 

215  re. 

220  if  hignrow  «  rowl  then  2*0 

225  re  rr 

25C  gosub  700  :  rei. 

2*0  gosub  7b0  :  re. 

250  gosub  820  :  re* 

255  reir, 

260  rowl  =  rowl  «  1  :  coll  '  coll  «  1  :  high  row  =  rowl  :  reir. increment  counte'S 

26 


re.  check  if  uniaue  solution 

re.  prints  out  srrjr 

re.  checks  for  row  with  highest  first  non-zero  entri 


--  eichenge  rows  if  necessary 

•-  sets  first  nor,  zero  entry  i  r,  rowl  to  1.0 

--  sets  first  entry  tc  zero  in  following  rows 


re. 
270  if  rowl 
280  gosub  7 
290  gosub  8 
295  re. 

500  r  e  »  

505  re. 
510  illcols 
520  for  i  « 
JJ0  for  j  = 
li)  * 
i 


()  rows  then  200  :      re. checks  if  lest  row  hes  beer,  reached 

60  :  re.  sets  first  non  zero  entry  in  final  row  to  1.0 

90  :  re.  final  printout  of  array 


calculates  the  unknown  i  values: 


540 


1*1  :  a  '    rows ! 

a  to  1  it  ep  - 1 

i  to  cols!  -  1 

illi)  -  ill  J*l )  *  All  i , j«l ) 


550  neit 
560  neit  1 
565  re. 
570  re.  --- 
575  re. 
580  for  k  « 
590  print 
400  neit  k 
410  end 
420  re.  --- 
425  re. 

450  re.  

4J5  re. 
500  for  i  « 
510  for  J  « 
520  if  i  * 
550  for  k  « 
540  if  Alii 
550  nilk )* 
560  neit  k 


prints  out  the  i  values: 


1  to  rows! 
'  k-1  "  » 


il  i  k  ) 


end  of  «ain  program 
begin  subroutines  -- 


1  to  rows ! 

1  to  rows ! 
j  then  610 

1  to  cols ! 
,k  )  »  0  then  560 

Al  I  J,  k  )  /Al  I  1  ,k  ] 
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570    for    ■    «    1    to    colt!    -    1 

580    if    nil*!     ()    nil.o]     then    610 

590    ntit    • 

600  print  "EQUATIONS  HOT  INDEPENDENT  ...  NO  UN10UE  SOLUTION  ...  ONE  SOLUTION  F0LL0US' 

610  next  J 

620  ntit  1 

630  return 

635  re* 

640  re»  subroutine  to  check  for  row  with  highest  first  non-zero  entry: 

6S0  for  1  »  rowl  to  row*!  -  1 

660  if  ■bslAKi.coll  )  )  )  absf  All  1»  I  .col  1  )  )  then  680 

670  highrow  ■  1  ♦  1 

680  neit  i 

690  return 

695  re»  »unroutine  to  eichenge  row*  and  put  highest  first  entry  first: 

700  for  1  >  coll  to  cols! 

710  Bt(i)  •  AKrowl,  1) 

720  Allrowl.i)  •  AKhlghrow,  i) 

730  Allhlghrow.  i)  •  Bl I  I  I 

740  neit  1 

750  return 

755  re*  subroutine  to  set  first  entry  in  first  row  to  1.0: 

760  Rl  '  AKrowl  .colli 

770  for  k  «  col  1  to  col s ! 

780  if  AKrowl, coll)  *  0  then  310 

79C  AKrowl. k)  •  AKrowl,  k)/RI 

791  print  k 

600  neit  k 

810  return 

815  re*  subroutine  to  set  other  first  entries  to  lerc: 

B2C  for  j  -    rowl  »  1  to  rows! 

830  PI  '    All j.col  1  l/AI(rowl,col 1 ) 

8iC  for  I  r  coll  to  cols! 

850  Allj.H  *  Al(j,k)  -  AKrowl, kl'PI 

860  neit  k 

870  ne»  t  J 

880  return 

885  re*  subroutine  to  print  out  the  tilrii  A: 

890  for  j  •  1  to  rows! 
900  f  or  k  *  1  to  col s ! 
910  print  All  J.k), 
920  neit  k 
930  print 
940  neit  i 

950  print  rowl;  ' - - - - 

9fc0  return 

970  print  "NOT  ENOUGH  E0UATI0NS  ...  HANT  SOLUTIONS  ARE  POSSIBLE' 

980  return 

990  end 
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PROGFAn    AZn 

*    -    6378506    4 

I    =    6356583.8 

0PEN(UNIT*1 ,F1LE= 'a: in. T" .FORM*' FORMATTED' . 
♦ACCESSs'SEOUENTlAL1 1 

REUIND    1 

srjmuti"     is    fro»    point     1     1 0    Doint     2 

READ! 1.2)    ALAT2.AL0NJ 

READfl.21  ALATI.AL0N1 

FLAG  =  0.0 

CALL  INVfA1B.ALATl.ALONl.ALA'2.ALON2.FLAT,DIST.AZI 

AZ  =  A2  •  3.  U1592654 

AZ  s  kill.  U15926S4M8C.  0 

IF  IAZ-J6C.0)  10.10.20 
20  AZ  =  AZ  -  JtO.O 
10  URITEI',3)  DIST.A2 

»  =  ALA72/J. U1592t54'180. 0 

Y  =  ALON2/3.  U1592654'  160.  0 

V  =  A  L  A  T 1 / J  14 1 5 9 : < 54 ' 1 E  C  0 
U  =  ALON1 /J. 14) 552154* 180.0 
UR1TE  < '. 4  )  >.  Y 

UP  I  TEC.  51  V.L 
2  FORMAT  fF9.e.TJ2.F9.el 


J  FOrr.n  1 

■DISTANCE 

=  "  .  F  J  2.  J, 

IX.  "METERS". JX. ' 

'  A  2  I M  'J  T  H  = 

".Fir 

.  '  I  N  0  F  T  H 

I' ) 

4  FOFMtTI 

'  P  C  J  K  T  2  :  ' 

. 2X, *LAT : ' 

.  2  >  .  F  E  : .  "  N  '  .  J  t 

'LPN  :  '  .  2> 

,  F  9 .  5 . 

5  F  OP'-*!  1 

'  P  0 1  N  T  1  :  ' 

.  ?)  ,  "LAT  ;  • 

.  2* . F£  5.  "N" , 3». 

"L0\:".2> 

,r;  5 . 

ENO 

SUBROUTINE  INV  [  EOF  AC'.  PORAP,  RLT  1  .  RLN1  .  RL  T2  .  RLN2.  CFLC-.  t'l  5T  .  A2  1 

i r  i  c  f  l  c- 1  i  o ; .  i  o  i  .  i  o ; 
ioi  fla'  =  i  g-forac  ec-fat 

FA.M:  -     FLA'*FLA" 
F  ;  =  F  A  L  A  T  2  "  1  .  2 I 

f;  =  falatj'o.j 

F3  -     FALA72'C.:r 
F  4  =  FALAT2'D.1IJ 

F  E  =  F  A  L  A  7  2  '  0  .  0 1 ."  E 

Ft   =  FtLAl;.FLA' 

r~    -    F  6  *  1  .  0 

F  E  =  F  6  •  0 .  5 

PI  =  3. 141592654 

TU0F1  =  b. 28318530' 

C  F  LC-  =  1.0 
10:  EETA]  =  ATANt ( 1 . O-FLAT 1  *  SI N( RL T 1 )/COS(RLTl)l 

5EETA1  =  S 1 N I e  E  T  A 1 ) 

CEETAl  =  COSIEETAl  1 

EETA:  =  ATANI  (  1  .C-FLAT  ) ' S IN ( RLT 2 ) /COS ( RLT 2 ] ) 

SEETA2  =  5 1 N ( B  E  T  A  2 1 

C  EE  T  A2  :  C  0  S I P  E  T  A  2 ) 

DELL  =  RLN1-FLN2 

ADELL  =  AESIOELL) 

IF  I  ADELL-FI  I  lOt  .  ICE  .  ICE 
105  ACELL  =  TUOFI-ACELL 
ICt  SIDEL  =  SIN t  ADELL 1 

CODEL  -    COSlADELLl 

A  =  SEE1A1*SEETA2 

e  =  cbetai'ceeta; 

C0FH1     =    A.fCODEL 
SIFHI     = 
650F1[(SlDEL,CEE1A:)-'2'l?EETA2"CErTAl-J.eET*l*CBETA:*CC'C'ELl"": 
C    =    B'SIDEL/SIPH] 
En    =     1 .O-C'C 
P  H  ]     =     A  S  I  N  I  S  1  P  H  I  1 
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IF     (COPHI)    207,208.208 

207  PHI    ■    PI-PHI 

208  PHISO    «    PHI'PMI 
CSPH1    «    J.0/S1PHI 
CTPHI    «    C0PHI/5IPHI 
PSYCO    «    SIPHI'COPHI 
TERni     *    F7'PH] 

TERH2    «    *'(Ft,SIPHI-F2'PHIS0*CSPHI  ) 

TERHJ    •    £n,(F2,PHI50*CTPHl-FB*(PMNPSYC0)) 

TERfU    »    A'A'F2,PSYC0 

TERH5    «    Ef1,En,(F5,(PHNPSYCO)-F2,PMIS0,CTPMl-F*,PSYCO,C0PHI,COPHII 

TERH6    «    »'Ef1*F2'  IPH1SO'CSPH1«PSYCO'COPH1  ) 

Dl  ST    «    PORAD' (TERHHTERn2«TERnJ-TERru«TERn5<TERH6) 

TERf17    «    Fb'PHI 

TERMS  ■  A,(F2,SIPHI«FALAT2*PH1S0,CSPHI) 

TERM9  •  EH'(FJ'F3'PSYC0«FALAT2,PH1S0,CTPHI-F1'PHI) 

2LAH  •  C*(TERM7-TERN8«TERH9)«ADELL 

CTAZ=(SBETA2'CBETA1-C0S(ZLAHZ),SBETA1,CBETA2)/(SIN(ZLAH),CBETA2) 

IF  (CTA2)  210,209.210 

209  CTAZ    «    0.00000005 

210  A2    *    ATAN( 1.0/CTA2) 
IF     (DELL)     21S.214.2U 

2U     IF     IDELL-PI  )    211.212.212 

211  IF     (CTAZ)    220, 221 . 221 
220    AZ    *    AZ«P] 

GOTO    221 

212  IF  (CTAZ)  217,218,218 
215  IF  (D£LL«PI )  211,211.21b 


21t 

IF  ICTAZ)  217,218.2 

217 

AZ  -  Pl-AZ 

GOTO  221 

218 

AZ  =  TUOPI-AZ 

221 

AZ  •  AZ.FI 

AZ  =  AZ-TUOPI 

If  (AZ)  213, 219. 219 

21  J 

AZ  =  AZ«TUOFI 

219 

RETURN 

END 
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APPENDIX  D 


122°00' 


121°50' 


33 


INITIAL  DISTRIBUTION  LIST 


1.  Defense  Technical  Information  Center 
Cameron  Station 

Alexandria,  VA  22314 

2.  Library,  Code  0142 
Naval  Postgraduade  School 
Monterey,  CA  93943 

3.  Chairman  (Code  68Co) 
Department  of  Oceanography 
Naval  Postgraduade  School 
Monterey,  CA  93943 

4.  Prof.  Stevens  P.  Tucker 

Depertment  of  Oceonography  (Code  68Tx) 
Naval  Postgraduade  School 
Monterey,  CA  93943 

5.  Director  Naval  Oceanography  Division 
Naval  Observatory 

34th  and  Massachusetts  Avenue  NW 
Washington,  DC  20390 

6 .  Commander 

Naval  Oceanography  Command 

NSTL  Station 

Bay  St.  Louis,  MS  39522 

7.  Commanding  Officer 

Naval  Oceanographic  Office 

NSTL  Station 

Bay  St.  Louis,  MS  39522 

8.  Commanding  Officer 

Fleet  Numerical  Oceanography  Center 
Monterey,  CA  93943 

9.  Commanding  Officer 

Naval  Ocean  Research  and  Development 

Activity 

NSTL  Station 

Bay  St.  Louis,  MS  39522 

10.   Commanding  Officer 

Naval  Environmental  Prediction 
Research  Facility 
Monterey,  CA  93943 
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11.  Chairman,  Oceanography  Department 
U.S.  Naval  Academy 

Annapolis,  MD  21402 

12.  Chief  of  Naval  Research 
800  N.  Quincy  Street 
Arlington,  VA  22217 

13.  Director  (Code  PPH) 
Defense  Mapping  Agency 

Bldg.  56,  U.S.  Naval  Observatory 
Washington,  DC  20305 

14.  Director  (Code  HO) 

Defense  Mapping  Agency  Hydrographic 
Topographic  Center 
6500  Brookes  Lane 
Washington,  DC  20315 

15.  Director  (Code  PSD-MC) 
Defense  Mapping  School 
Ft.  Belvoir,  VA  22060 

16.  Director,  Charting  and  Geodetic 
Services  (N/CG) 

National  Ocean  and  Atmospheric 
Administrat  ion 
Rockville,  MD  20852 

17.  Chief,  Program  Planning,  Liaison 
and  Training  (NC2) 

National  Oceanic  and  Atmospheric 

Administration 

Rockville,  MD  20852 

18.  Chief,  Nautical  Charting  Division  (N/CG2) 
National  Oceanic  and  Atmospheric 
Administration 

Rockville,  MD  20852 

19.  Chief,  Hydrographic  Surveys  Branch 
(N/CG24 ) 

National  Oceanic  and  Atmospheric 

Administration 

Rockville,  MD  20852 

20.  Director,  Pacific  Marine  Center 
(N/MOP) 

National  Ocean  Service,  NOAA 
1801  Fairview  Avenue  East 
Seattle,  WA  98102 


85 


21.  Director,  Atlantic  Marine  Center 
(N/MOA) 

National  Ocean  Service,  NOAA 
439  W.  York  Street 
Norfolk,  VA  23510 

22.  Cdr.  Kurt  J.  Schnebele,  NOAA 
Dept.  of  Oceanography  (Code  68SN) 
Naval  Postgraduate  School 
Monterey,  CA  93943 

23.  Capt.  Michael  Iacona  (USNR  ret) 
2304  Roslyn  Ave. 

District  Hts . ,  MD  20747 

24 .  Comandant 
G-NRN-1 

2100  SW  2nd  St.  SW 
Washington,  DC  20593-0001 

25.  Defense  Mapping  Agency 
Hydrographic/Topographic  Center 
Code  HCNN,  Attn:  L.  Phillips 
6500  Brookes  Lane 
Washington,  DC  20315 

26.  IHO/FIG  International  Advisory  Board 
International  Hydrographic  Bureau 
Avenue  President  J.F.  Kennedy 
Monte  Carlo,  Monaco 

27.  Hydrographer  of  the  Navy 
Ministry  of  Defense 
Hydrographic  Department 
Taunton,  Somerset 
TAI29N  England 

28.  Hydrographer,  Royal  Australian  Navy 
Hydrographic  Office 

Box  1332 

North  Sydney  2060 

New  South  Wales,  Australia 

29.  Dominion  Hydrographer 
Canadian  Hydrographic  Service 
615  Booth  Street 

Ottawa,  Ontario 
Canada  K1AOE6 

30.  Kapala  Jawatan  Hidro  Oseanografi 
Jalan  Gunung  Sahari  -  87 
Jakarta,  Indonesia 
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31.  Direccion  General  De  Oceanografia 
Departamento  De  Hidrograf icia 
Ave.  COYOACAN  131 

Mexico  12,  D.F. 

32.  Director 

Instituto  Hidrografico 
Rua  Das  Trinas,  49 
Lisboa,  Portugal 

33.  Civil  Engineering  Department 

Ports  and  Lighthouses  Administration 
Ras  Eltin  -  Gate  1 
Alexandria,  Egypt 

34.  Hydrographic  Service 
Athens  BST  902,  Greece 

35.  Chief,  Hydrographic  Surveys  Branch 
Atlantic  Marine  Center,  NOAA 

439  West  York  St. 
Norfolk,  VA  23510 

36.  Commanding  Officer 

Naval  Western  Oceanography  Center 

Box  113 

Pearl  Harbor,  HI  96860 

37.  Commanding  Officer 

Naval  Oceanography  Command  Center,  Rota 

Box  31 

FPO  San  Francisco,  CA   96540 

38.  Library 

Moss  Landing  Marine  Lab 
California  State  Colleges 
Sandholdt  Road 
Moss  Landing,  CA  95039 

39.  Library 

School  of  Oceanography 
Oregon  State  University 
Corvallis,  OR  97331 

40.  Library 

Dept.  of  Oceanography 
University  of  Washington 
Seattle,  WA   98105 

41.  ARNAV  Systems  Inc.  (attn:  Mr.  Walter  Dean) 
4740  Ridge  Dr.  NE 

Salem,  OR  97303-0012 
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42.   Regional  Director 

U.S.  Dept .  of  the  Interior 
Minerals  Management  Service 
1340  W.  6th  St. 
Los  Angeles,  CA  90017 
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